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We review the basic mechanisms of neutrino mass generation and the corresponding structure 
of the lepton mixing matrix. We summarize the status of three-neutrino osciUation parameters 
as determined from current observations, using state-of-the-art solar and atmospheric neutrino 
fluxes, as well as latest experimental data as of September 2007. We also comment on recent 
attempts to account for these results and to understand flavour from first principles. We discuss 
extensively the prospects for probing the strength of CP violation in two near term accelerator 
neutrino oscillation experiments, T2K and NOz^A, as well as possible extensions such as T2KK and 
a second large off-axis detector near the NOz^A detector. We also briefly discuss the possibility 
of probing the effect of Majorana phases in future neutrinoless double beta decay searches and 
discuss other implications of leptonic CP violation such as leptogenesis. Finally we comment on 
the issue of robustness of the current oscillation interpretation and possible ways of probing for 
non-standard neutrino interactions in precision oscillation studies. 
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1 Introduction 

CP violation has been discussed phenomenologically since the sixties. Gauge theories naturally account 
for its presence, although leave its ultimate origin unanswered [1]. Understanding the origin of CP 
violation from first principles constitutes one of the central challenges of theoretical elementary particle 
physics. With the historic discovery of neutrino oscillations [2, 3, 4, 5, 6] the issue of leptonic CP 
violation has also come to the center of the agenda of the particle and nuclear physics communities. 
The existence of CP violation in the lepton sector is expected in gauge theories of neutrino mass. The 
main difference with respect to CP violation in the quark mixing matrix is the appearance of new phases 
associated to the Majorana nature of neutrinos and/or with the admixture of SU{3) SU{2) U{1) 
singlet leptons in the charged current weak interaction [7, 8]. The latter also leads to the effective 
violation of unitarity in the lepton mixing matrix describing neutrino oscillations. Irrespective of what 
the underlying origin of neutrino mass may turn out to be, Majorana phases are a generic feature of 
gauge theories that account for the smallness of neutrino mass through the feebleness of lepton (or B-L) 
number violation. This includes two alternative classes of gauge theories, that differ by the scale at 
which L is broken and neutrino masses arise. In Sec. 3 we give a very sketchy summary of the various 
ways to endow neutrinos with mass, and in Sec. 4 we describe the basic structure of the lepton mixing 
matrix that follows from theory. 

The analysis of current neutrino oscillation experiments is given in Sec. 5 within the simplest CP- 
conserving three-neutrino mixing pattern, leaving aside the LSND and Mini-Boone data [9, 10, 11]. 
We summarise the status of neutrino mass and mixing parameters [12] as determined from current 
neutrino oscillation data [2, 3, 4, 5, 6] [14, 15]. In addition to a determination of the solar angle ^12, the 
atmospheric angle 623 and the corresponding mass squared splittings, the data place a constraint on the 
last angle in the three-neutrino leptonic mixing matrix, ^13. Necessary inputs of such interpretation are 
the solar and atmospheric neutrino fluxes [16, 17], the neutrino cross sections and response functions, as 
well as the accurate description of neutrino propagation in the Sun and the Earth, taking into account 
matter effects [18, 19]. 

It is well known that CP violation should manifest itself in neutrino oscillation experiments [20]. 
In Sec. 6 we review the basic theoretical features of CP violation in neutrino oscillation probabilities. 
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while in Sec. 7 we proceed to a summary of the prospects for probing ^13 and CP violation in neutrino 
oscillations. We focus on the case of accelerator neutrino experiments and give a detailed discussion of 
the potential of the T2K and NOi/A experiments and their possible extensions such as T2KK and a 
second large off-axis detector near the NOi/A detector. 

Lepton number violating processes such as neutrinoless double beta decay [21, 22] are briefly dis- 
cussed in Sec. 8. In particular, searching for (5(5^^ constitutes a very important milestone for the 
future, as this will probe the fundamental nature of neutrinos, irrespective of the nature of the neutrino 
mass-generation mechanism [23]. Within the simplest mass mechanism the search for /3/3oi/ probes not 
only the absolute scale of neutrino mass but is also sensitive to CP violation induced by the so-called 
Majorana phases [7], inaccessible in conventional oscillations [24, 25, 26, 27]. 

In Sec. 9 we also briefly comment on the robustness of the current oscillation interpretation and 
point out the interest in probing for non-standard neutrino interactions in future precision oscillation 
studies. 



2 Kinematics of neutrino mass 
2.1 Majorana and Dirac masses 

Massive fermions can either be Dirac or Majorana. If they carry electric charge, wc have no choice, they 
must be Dirac. Electrically neutral fermions, like neutrinos (or supersymmetric " inos"), arc expected 
to be Majorana-type on general grounds, irrespective of how they acquire their mass (see Sec. 3). 
Phcnomcnological differences between Dirac and Majorana neutrinos arc tiny for most processes, such 
as neutrino oscillations: first because neutrinos arc known to be light and, second, because the weak 
interaction is chiral, well described by the V-A form. Nevertheless it is basic, in particular, for the issue 
of CP violation and, for this reason, it will be reviewed here. 

The most basic spin 1/2 fermion corresponding to the lowest representation of the Lorentz group is 
given in terms of a 2-component spinor p, with the following free Lagrangean [7] 

Tfi 

J-'M = -ip^c^nd^p - —p^(T2p + H.C. (1) 
where we use the 2 x 2 cr matrices, with (jj being the usual Pauli matrices and (74 = —i I, I being 

— * 

the identity matrix, in Pauli's metric conventions, where a.h = a^h^j^ = a ■ b + 0464, 04 = iag. Under a 
Lorentz transformation, x — > Ax, the spinor p transforms as p ^ S{A)p{A^^x) where S obeys 

S^a^^S = A^^a, . (2) 

The kinetic term in Eq. (1) is clearly invariant, and so is the mass term, as a result of unimodular 
property det S — 1. However, the mass term is not invariant under a phase transformation 

P - e'-p. (3) 
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The equation of motion following from Eq. (1) is 



-icF^d^p = ma2p*. (4) 

As a result of the conjugation and Clifford properties of the cr-matrices, one can verify that each 
component of the spinor p obeys the Klein-Gordon wave-equation. 

In order to display clearly the relationship between our theory, Eq. (1) and the usual theory of a 
massive spin 1/2 Dirac fermion, defined by the familiar Lagrangean 

Cv = - m (5) 

where by convenience we use the chiral representation of the Dirac algebra ^i/y^ + Ivlfx — 2 5^1, in which 
75 is diagonal, 

In this representation the charge conjugation matrix C obeying 

= -C (7) 
= C-^ (8) 
l,C^-^l (9) 

(10) 

is simply given in terms of the basic conjugation matrix a"2 as 



A Dirac spinor can then be written as 

so that the corresponding charge- conjugate spinor ^'^ = is the same as ^'/j but exchanging and 

X, i- e. 

*i,=f ^). (13) 

A 4-component spinor is said to be Majorana or self-conjugate if = which amounts to setting 
X — 4>- Using Eq. (12) we can rewrite Eq. (5) as 

2 2 

Cv^-iY^ Pi(^ndf,pa - y XI ^«^2Pa + H.C. (14) 

a=l a=l 



where 



1 



(P2 + 



X 



^/2 



1 



{p2-ipi), 



(15) 



V2 



are the left handed components of and of the charge-conjugate field ^'^j respectively. In this way the 
Dirac fermion is shown to be equivalent to two Majorana fermions of equal mass. The U{1) symmetry 
of the theory described by Eq. (5) under ^ e'"^i5 corresponds to continuous rotation symmetry 
between pi and p2 



which result from the mass degeneracy between the p's, showing that, indeed, the concept of fermion 
number is not basic. 

2.2 Quantization 

The mass term in Eq. (1) vanishes unless p and p* are anti-commuting, so wc consider the Majorana 
fermion, right from the start, as a quantized field. The solutions to Eq. (1) can easily be obtained in 
terms of those of Eq. (5), which are well known. The answer is 



operators obey canonical anti-commutation rules and, like the m's and f 's, depend on the momentum 
k and helicity label r. The expression in Eq. (16) describes the basic Fourier expansion of a massive 
Majorana fermion. It differs from the usual Fourier expansion for the Dirac spinor in Eq. (17) in two 
ways, 

• the spinor is two-component, as there is a chiral projection acting of the it's and v^s 

• there is only one Fock space, particle and anti-particle coincide, showing that a massive Majorana 
fermion corresponds to one half of a conventional massive Dirac fermion. 

The It's and v's are the same wave functions that appear in the Fourier decomposition the Dirac field 



Pi cos 6pi + sin 6p2 
P2 — — sin 9pi + cos 9p2 




^^ = (27r)-^/2 / d=^A;^(^)V2[e^'=-X(A;)it,(A;) + e-^*^-^6t(A;)i;,(A;)] 



(17) 



r=l 

Using the helicity eigenstate wave-functions. 



a ■ kuf{k) = ± I A; I u^{k), 
a-kvt{k) = T\ k\ v^{k), 



(18) 
(19) 
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one can show that, out of the 4 hnearly independent wave functions uf^{k) and v^{k), only two survive 
as the mass approaches zero, namely, ul{k) and v'^{k) [28]. This way the Lee- Yang two-component 
massless neutrino theory is recovered as the massless limit of the Majorana theory. 
Two independent propagators follow from Eq. (1), 

< I p{x) p*{y) I >= ia^d^I\F{x -y;m), (20) 
< I p{x) p{y) I >= m (72 Af(x - y; m), (21) 

where /S.f{x — y;m) is the usual Feynman function. The first one is the "normal" propagator that 
intervenes in total lepton number conserving (AL = 0) processes, while the one in Eq. (21) describes 
the virtual propagation of Majorana neutrinos in AL = 2 processes such as neutrinoless double-beta 
decay. 

2.3 CP properties 

The Lagrangean in Eq. (1) can be easily generalized for a system of an arbitrary number of Majorana 
neutrinos, giving 

n 1 " 

CM = -tYl Pi^M^MPa - 2 E Ma0l(T2Pp + H.C. (22) 

a=l a,/3=l 

where the sum runs over the "neutrino-type" indices a and j3. By Fermi statistics the mass coefficients 
Map must form a symmetric matrix, in general complex. This matrix can always be diagonalized by a 
complex n X n unitary matrix U [7] 

Mdiag = U^MU . (23) 

When M is real (CP conserving) its diagonalizing matrix U may be chosen to be orthogonal and, in 
general, the mass eigenvalues can have different signs. These may be assembled as a signature matrix 

r)^diag{+, +,..., -,-,..) (24) 

For two neutrino types there are two classes of models, one with rj — diag{+, — ) and another char- 
acterized by = diag{+, +). The class with rj — diag{+, — ) contains as a limit the case where the 
two fermions make up a Dirac neutrino. Depending on whether the two neutrinos are active-sterile or 
active-active the limit is called quasi- [29, 30] or pseudo-Dirac neutrino [31]. Oscillations between the 
quasi-Dirac neutrino components potentially change the predictions of Big Bang Nucleosynthesis while 
oscillations between the components of a pseudo-Dirac neutrino do not. 

Note that one can always make all masses positive by introducing appropriate phase factors in the 
wave functions, such as the factors of i in Eq. (15). When interactions are added (see Sec. 4) these signs 
become physical. As emphasized by Wolfenstein, these signs and the corresponding CP phases play an 
important role in the discussion of neutrinoless double beta decay [32] . 
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SU{3)0 SU{2)0U{1) 


f — (u 1 V 


(1,2,-1) 




(1,1,2) 


Qa (^a, da) 


(3,2,1/3) 


K 


(3, 1, -4/3) 


d" 


(3,1,2/3) 


$ 


(1,2,1) 



Table 1: Matter and scalar multiplets of the Standard Model 

3 The origin of neutrino mass 

The fifteen basic building blocks of matter listed in Table 1 are all 2-component sequential "left-handed" 
chiral fermions, one set for each generation. Parity violation in the weak interaction is accounted 
for "effectively" by having "left" and "right" fermions which behave differently with respect to the 
SU (3) ® SU (2) f/(l) gauge group. In contrast to charged fermions, neutrinos come only in one chiral 
species. It has been long noted by Weinberg [33] that one can add to the Standard SU {3)^SU (2) (8)f/(l) 
Model (SM) an effective dimension-five operator O = A££<l>$, where i denotes a lepton doublet for 
each generation and $ is the SM scalar doublet. 




Figure 1: Dimension five operator responsible for neutrino mass. 

Once the electroweak symmetry breaks through the nonzero vacuum expectation value (vev) ($) 
Majorana neutrino masses oc ($)^ are induced, in contrast to the masses of the charged fermions which 
are linear in ($). This constitutes the most basic definition of neutrino mass, in which its smallness 
relative to the masses of the SM charged fermions is ascribed to the fact that O violates lepton number 
by two units (AL = 2) whereas the other fermion masses do not. Note that this argument is totally 
general and holds irrespective of the underlying origin of neutrino mass. From such general point of view 
the emergence of Dirac neutrinos would be a surprise, an "accident" , justified only in the presence of 
a fundamental lepton number symmetry, in general absent. For example, neutrinos could naturally get 
very small Dirac masses via mixing with a bulk fermion in models involving extra dimensions [34, 35, 36]. 
Barring such very special circumstances, gauge theories expect neutrinos to be Majorana. 

Little more can be said from first principles about the mechanism giving rise to the operator in 
Fig. 1, its associated mass scale or its flavour structure. For example, the strength A of the operator 
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O may be suppressed by a large scale Mx in the denominator (top-down) scenario, leading to 



where Aq is some unknown dimensionless constant. Gravity, which in a sense "belongs" to the SM, 
could induce the dimension-five operator O , providing the first example of a top-down scenario with 
Mx — Mp, the Planck scale. In this case the magnitude of the resulting Majorana neutrino masses are 
too small to be relevant in current searches. 

Alternatively, the strength A of the operator O may be suppressed by small parameters (e.g. scales, 
Yukawa couplings) in the numerator and/or loop- factors (bottom- up scenario). Both classes of scenarios 
are viable and allow for many natural realizations. While models of the top-down type are closer to the 
idea of unification, bottom-up schemes are closer to experimental verification. 

Models of neutrino mass may also be classified according to whether or not additional neutral heavy 
states are present, in addition to the three isodoublet neutrinos. As an example, such leptons could be 
SU{3) ® SU{2) ® U{1) singlet "right-handed" neutrinos. In what follows we classify models according 
to the mass scale at which O is induced, namely bottom-up and top-down scenarios. 

3.1 Seesaw- type neutrino masses 

The most popular top-down scenario is the seesaw. The idea is to generate the operator O by the 
exchange of heavy states. The smallness of its strength is understood by ascribing it to the violation of 
lepton number at a high mass scale, namely the scale at which the heavy states acquire masses. 

3.1.1 The majoron seesaw 

The simplest possibility for the seesaw is to have ungauged lepton number [8]. It is also the most 

general, as it can be studied in the framework of just the iS'f/(3) ® SU{2) ® U{1) gauge group. Such 
"1-2-3" scheme is characterized by 5*^7(3) ® SU{2) ® U{1) singlet, doublet and triplet mass terms, 
described by the matrix [7, 8] 



in the basis z/l, '^l, corresponding to the three "left" and three "right" neutrinos, respectively. Note that, 
though symmetric, by the Pauli principle, M.^ is complex, so that its Yukawa coupling sub-matrices 
Fj, as well as Y3 and Yi are complex matrices denoting the relevant Yukawa couplings, the last two 
symmetric. 

Such SU{?>) ® SU{2) ® U{1) seesaw contains singlet, doublet and triplet scalar multiplets, obeying 
a simple "1-2-3" vev seesaw relation of the type 



This follows simply from the minimization condition of the S'f/(3) ® SU{2) U{1) invariant scalar 
potential and arises in a wide variety of seesaw type models, as reviewed in [37, 38]. It implies that the 




(25) 



V3V1 ~ with ■Ui » ■U2 » V3 



(26) 
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triplet vev f 3 — as the singlet vev Vi grows. Here V2 = ($) denotes the SM Higgs doublet vev. Small 
neutrino masses are induced either by heavy SU{3) ® SU{2) ® f/(l) singlet "right-handed" neutrino 
exchange (type I) or the smallness of the induced triplet vev that follows from heavy scalar exchange 
(type II), as illustrated in Fig. 2. The matrix Ai,^ is diagonalized by a unitary mixing matrix Uiy, 




Figure 2: Two types of seesaw mechanism. Left: right-handed neutrino exchange (type-I), right: heavy 
SU{3) ® SU{2) ® U{1) triplet exchange (type-II). 



UjMM. = dmg{m,,M,), (27) 

yielding 6 mass eigenstates, including the three light neutrinos with masses rrii, and three two-component 
leptons. The light neutrino mass states z/j are given in terms of the flavour eigenstates via the unitary 
matrix 11^, [7, 8] 

6 

l^i = J^iUuhaUa. (28) 

a=l 

The effective light neutrino mass, obtained this way is of the form 

~ ^3^3 - Y,Y,'V^ (29) 

Vl 

The diagonalization matrices can be worked out explicitly as a perturbation series, see Ref. [8]. 

There is an important feature of such "1-2-3" seesaw namely, that since lepton number is ungauged, 
there is a physical Goldstone boson associated with its spontaneous breakdown, the Majoron [39]. It is 
often argued that, due to quantum gravity effects the associated majoron will pick up a mass. It has 
been shown that, a massive majoron with mass in the fraction of keV range can provide the observed 
dark matter of the Universe [40]. Such majoron decaying dark matter offers an alternative to the 
supersymmetric dark matter scenario and may lead to interesting phenomenological implications. 

3.1.2 Left-right symmetric seesaw 

A more elegant setting for the seesaw is a gauge theory containing B-L as a generator, such as SU (3) ® 
SU{2)l ® SU{2)n ® U{1)b^l or the unified models based on ^O(IO) or E{6) [41, 37, 42]. For example 
in 50(10) each matter generation is naturally assigned to a 16 (spinorial in 50(10)) so that there are 



10 



16 . 16 . 10 as well as 16 . 126 . 16 terms generating Dirac and Majorana neutrino mass terms, 
respectively, leading to the neutrino mass matrix 



Yl (Al) ($) 

Yr{Ar) 



M.-\ i m 



in the basis vl, i/£, corresponding to the "left" and "right" neutrinos, respectively, where Yl and Yr 
denote the Yukawas of the 126 of 5'O(10), whose vevs {Al^r) give rise to the Majorana terms, while 
Yl, is the Dirac Yukawa coupling in 16 . 16 . 10. The matrix Y^, is an arbitrary complex matrix 
in flavour space, while Yl and Yr are complex symmetric 3x3 matrices that correspond to Yi and 
1^3 of the simplest SU{3) (8) SU{2) (8) U{1) model. Small neutrino masses are induced either by heavy 
SU {S)<SiSU (2) (8)C/(1) singlet "right-handed" neutrino exchange (type I) or heavy scalar boson exchange 
(type II) as illustrated in Fig. 2. The matrix A4i, is diagonalized by a unitary mixing matrix Ui, as 
before. The diagonalization matrices can be worked out explicitly as a perturbation series, using the 
same method of Ref. [8]. This means that the exphcit formulas for the 6x6 unitary diagonalizing 
matrix U given explicitly in Ref. [8] also hold in the left-right case, provided one takes into account 
that vi — > (Ar) and vs — > (A^). The effective hght neutrino mass, obtained this way is of the form 

Yl (Al) - Y^Yr'^Y/^. (31) 

We have the new vev seesaw relation 

(Al) (An) {^f , (32) 

which naturally follows from minimization of the left-right symmetric scalar potential, together with 
the vev hierarchy [37] 

{Al) < ($) < (A,^) . (33) 

This imphes that both type I and type II contributions vanish as (Ar) oo. The structure of the 
seesaw is exactly the same as before, enabling one to employ the same perturbative diagonalization 
method in [8]. 

The new insight that is obtained this way is that now one can arrange for the breakdown of parity 
invariance to be spontaneous, so that smallness of neutrino masses gets correlated to the observed maxi- 
mality of parity violation in low-energy weak interactions, as stressed by Mohapatra and Senjanovic [37]. 
However this is hardly relevant phenomenologically in view of the large value of the B-L scale needed 
both to fit the neutrino masses, as well as unify the gauge couplings. Another important difference with 
the SU (3) (8) SU (2) (8) U (1) seesaw case is the absence of the Majoron, now absorbed as the longitudinal 
mode of the gauge boson corresponding to the B-L generator. 

3.1.3 Double seesaw 

Nothing is sacred about the number of (anomaly-free) gauge singlet leptons Si in the SM [7] or of 
singlets outside the 16 in SO{10) or the 27 in E{6) [43]. New important features may emerge when the 
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seesaw is realized with non-minimal lepton content. Here we mention the seesaw scheme suggested in 
Ref. [43] with E{6) motivations [44]. The model extends minimally the particle content of the SM by 
the sequential addition of a pair of two-component SU (3) ® SU{2) ® U{1) singlet leptons, v^, Si, with i 
a generation index running over 1, 2, 3. In the z/, u"^, S basis, the 9x9 neutral leptons mass matrix Ai^, 
is given as 

/ rj(<i>) A 

Mu= ¥,{<!>) , (34) 

^0 M J 

in the basis i>i, z/£, Sl, where is an arbitrary 3x3 complex Yukawa matrix, M and fi are SU{3) ® 
SU{2) U{1) singlet complex mass matrices, fi being symmetric. Notice that it has zeros in the vl-^l 
and vl-vl entries, a feature of several string models [44]. 

For ^ M one has to first approximation that the Si decouple leaving the simpler seesaw at scales 
below that. In such a "double" seesaw scheme the three light neutrino masses are determined from 



(35) 



The mass generation is illustrated in Fig. 3. This formula can be readily derived from the method given 




Figure 3: "Double" and "inverse" seesaw mechanism. 



in Ref. [8] and decoupling the heavy states in two steps, first S then u'^ since /i 3> M. A new feature 
is that there are two independent scales, fi and M, with the B-L symmetry broken only by the largest 
scale /i. It is completely natural also to consider the case where /i instead of large is very small [43], 
even smaller than the electroweak vev. This low fi case will be considered in Sec. 3.2.3, and opens new 
phenomenological possibilities, discussed in Sec. 3.3. 

Irrespective of what sets its scale, the entry fi may be proportional to the vev of an SU{3) ® SU (2) ® 
f/(l) singlet scalar, in which case the model contains a singlet majoron. 



3.1.4 Unconventional seesaw 

There are many types of seesaw. More important than keeping track of the taxonomy of schemes (type 
I [41, 37, 42], type II [7, 8], type III [45, 46, 47], etc.) is understanding that the seesaw is not a theory 
but a mechanism that allows for many possible realizations. Schemes leading to the same pattern of 
neutrino masses may differ in many other respects. 
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As an example of extended seesaw models, let us consider one that has recently been suggested [48]. 
It belongs to the class of supersymmetric SO (10) models with broken D-parity. In addition to the three 
right-handed neutrinos of the standard seesaw model, it contains three sequential gauge singlets SiL 
with the following mass matrix 



( 

V iXL) 



F{xr) 

F^ (xb) ; 



(36) 



in the basis ul, i/£, Sl- Notice that it has zeros along the diagonal, specially in the vl-^l and vl-vl 
entries, thanks to the fact that there is no 126, a feature of several string-inspired models [44, 43]. The 
resulting neutrino mass is 



unif 



(37) 



where Munif is the unification scale, F and F denote independent combinations of Yukawa couplings of 
the SiL- One can see that the neutrino mass is suppressed by the unification scale Munif irrespective of 
how low is the B-L breaking scale. In contrast to all familiar seesaws, see e.g. Eq. (31), this new seesaw 
mechanism is linear in the Dirac Yukawa couplings Yy, as illustrated in Fig. 4. It is rather remarkable 



Xr 



5|C 



XL 



Xr 



1/ 



Figure 4: Unconventional seesaw mechanism. 

that one can indeed take the B-L scale as low as TeV without generating inconsistencies with gauge 
coupling unification [48]. 



3.2 Bottom-up neutrino masses 

There is a variety of models of neutrino mass where the operator O is induced from physics at accessible 
scales, TeV or less. The smallness of its strength is naturally achieved due to loop and Yukawa couplings 
suppression. Note also that lepton number violating parameters may appear in the numerator instead 
of the denominator of O with the smallness of its strength natural in t'Hofft's sense [49]. For example, 
in the inverse seesaw scheme of Sec. 3.2.3, one can consistently take /x to be small as the symmetry of 
the theory increases in the limit of vanishing /x, namely B-L is restored. 

3.2.1 Radiative models 

The first possibility is that neutrino masses are induced by calculable radiative corrections [50]. For 
example, they may arise at the two-loop level [51] as illustrated in Fig. 5. Up to a logarithmic factor 
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one has, schematically, 

2 /^\2 



^^-''[m-^) ^'ifF^"^* '''' 

in the limit where the doubly-charged scalar k is much heavier than the singly charged one. Here / 
denotes a charged lepton, / and h are their Yukawa coupling matrices and Yi denotes the SM Higgs 
Yukawa couplings to charged leptons and (a) is an SU{3) SU{2) U{1) singlet vev introduced in 
Ref . [52] . Clearly, even if the proportionality factor Aq is large, the neutrino mass can be made naturally 
small by the presence of a product of five small Yukawas and the appearance of the two-loop factor. A 
remarkable feature of the model is that, thanks to the anti-symmetry of the / Yukawa coupling matrix, 
one of the neutrinos is massless. 



Figure 5: Two- loop origin for neutrino mass. 



3.2.2 Supersymmetric neutrino masses 

Another very interesting alternative are models where low energy supersymmetry is the origin of neu- 
trino mass [53] . The intrinsically supersymmetric way to break lepton number is to break the so-called 
R parity. This could happen spontaneously, driven by a nonzero vev of an SU (3) ® SU (2) U{1) singlet 
sneutrino [54, 55, 56]. This way one is led to a very simple reference model which can be regarded as 
the minimal way to include neutrino masses into the MSSM. In this model R parity is violated only 
through an effective bilinear term [57]. Neutrino mass generation takes place in a hybrid scenario, with 
one scale generated at tree level by the mixing of neutralinos and neutrinos, induced by the sneutrino 
vevs, and the other induced by "calculable" radiative corrections [58]. Here the two blobs in each graph 

^3 

— 

Figure 6: Loop origin of solar mass scale. Atmospheric scale arises from a type-I weak-scale seesaw-like 
graph involving exchange of sleptons and squarks. 

denote AL = 1 insertions, while the crossed blob accounts for chirality flipping. The general form of 
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the expression is quite involved but the approximation 

- (i^) 

holds in some regions of parameters. Here A denotes the trilinear soft supersymmetry breaking coupling 
(not shown in the diagram), mo the scalar supersymmetry breaking mass, and Yd the suitable Yukawa 
coupling. The neutrino mass spectrum naturally follows a normal hierarchy, with the atmospheric scale 
generated at the tree level and the solar mass scale arising from calculable loops. 

3.2.3 Inverse seesaw 

Before closing this section we mention that there are also tree level neutrino mass schemes with naturally 
light neutrinos. One is the inverse seesaw scheme suggested in Ref. [44] from heterotic string motivations 
which led to the £'(6) gauge group. The model extends minimally the particle content of the SM by 
the sequential addition of a pair of two-component SU{^) ® SU{2) U{1) singlet leptons, uf, S^, with 
i denoting a generation index running over 1, 2, 3. 

In the v,v'^,S basis, the 9x9 neutral leptons mass matrix M.^, is the formally exactly the same as 
given in Eq. (34) [43]. Again the diagonalization follows the method in Ref. [8] but, in contrast to what 
is done for the large n regime, now one separates the full heavy sector consisting of three Quasi- Dirac 
states (six two-component leptons) at once. One obtains exactly the same light neutrino mass formula 
despite the important difference that now the entry /i is taken very small, e. g. <C Yy ($) <C M. 

As before Y^ and Af are arbitrary 3x3 complex Yukawa matrices, /x being symmetric due to the 
Pauli principle. Notice that for small /i neutrino masses vanish with yU, as we saw in Eq. (35), illustrated 
in Fig. 3. The fact that the neutrino mass vanishes as /i is just the opposite of the behaviour of the 
seesaw formulas in Eqs. (29) and (31) with respect to ^3 and (A^^), respectively; thus this is sometimes 
called inverse seesaw model of neutrino masses. The entry may be proportional to the vev of an 
SU{?)) ® SU{2) (g) U{1) singlet scalar, in which case spontaneous B-L violation leads to the existence of 
a majoron [59]. This would be characterized by a relatively low scale, so that the corresponding phase 
transition could take place after the electroweak transition [52] . The naturalness of the model stems from 
the fact that the limit when /i ^ increases the symmetry of the theory. One possible phenomenological 
implication would be the phenomenon of invisibly decaying Higgs boson [60, 61, 62]. In such schemes it 
will be crucial to take into account the existence of sizeable invisible Higgs boson decay channels in the 
analysis of experimental data on Higgs searches [63, 64]. Another possible implication is the existence 
of novel neutrino decay and annihilation processes that may be relevant in dense supernova media [65]. 
Other aspects of the phenomenology are mentioned in Sec. 3.3. 

3.3 Phenomenology of neutrino masses and mixings 

Clearly the first phenomenological implication of neutrino mass models is the phenomenon of neutrino 
oscillations, required to account for the current data. Before turning to the prospects for probing 
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neutrino oscillations with high precision at the future generation of long baseline oscillation experiments 
discussed in Sec. 6, let us comment on other manifestations of neutrino mass. 

First comes the issue of the absolute scale of neutrino masses. This will be tested directly in searches 
for tritium beta decay [66, 67] and cosmology [68, 69]. Some models do indeed suggest a lower bound on 
the absolute neutrino mass. For example the model in Ref. [70] gives m^, ^ 0.3 eV and therefore will be 
tested soon. Lastly, the nature (Dirac or Majorana) of neutrinos and the new sources of CP violation 
associated to its Majorana character can be scrutinized in neutrinoless double beta decay searches and 
other lepton number violation processes (see Sec. 8). What else is there? 

If neutrino masses arise a la seesaw this may be all. The simplest SO (10) seesaw has a drawback 
that the dynamics responsible for generating the small neutrino masses seems most likely untestable. In 
other words, beyond neutrino masses and oscillations, the model can not be probed phenomenologically 
at low energies, due to the large scale involved. 

However, in the presence of supersymmetry, the phenomenology of the seesaw mechanism can be 
richer. A generic feature of supersymmetric seesaw models is the existence lepton flavour violation 
decays like fi~ — > e~7, flavour violating tau decays (Fig. 7) and nuclear fi^ — conversion (Fig. 8). 
These can have accessible rates (Figs. 9 and 10) which depend not only on the seesaw mechanism, but 
also on the details of supersymmetry breaking and on a possible theory of flavour. The way these LFV 
processes arise through supersymmetry [71] is illustrated in Fig. 7. The existence of such loop effects 
leads to enhanced rates for flavour violating processes [72, 73, 74, 75, 76, 77] [78, 79]. We note that 
such LFV and/or CP violating effects may arise at the one-loop level from the exchange of relatively 
light neutral heavy leptons, even in the absence of supersymmetry [80, 81, 82, 83, 84]. 




Figure 7: Supersymmetric Feynman diagrams for I- — > Z^- 7. They involve the exchange of 
charginos(neutralinos) and sneutrinos (charged sleptons). 

As an illustration of the interplay of supersymmetry with heavy lepton exchange in engendering 
lepton flavour violation processes we consider the rates for the fi^ e~7 decay in the framework of 
the supersymmetric inverse seesaw model [78, 79]. Fig. 9 displays the dependence of the branching 
ratios for /i~ — e~ conversion in Ti (left) and fi^ e~7 (right) with the small neutrino mixing angle 
6'i3, for different values of 6*12 (black curve: 612 best fit value, blue bands denote 2a, 3cr, 4cr confidence 
intervals for the solar mixing angle 612). The inverse seesaw parameters are given by: M = 1 TeV and 
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/X = 30 eV. The light neutrino parameters used are from [12], except for ^13 which is varied as shown 
in the plots. The vertical coloured bands denote the exclusion areas from the 2(T, 3a and Aa limits 
on sin^6'i3, respectively. These lepton flavour violation rates may be testable in the new generation of 
upcoming experiments. For large M the estimates are similar to those of the standard supersymmetric 
seesaw. 




Figure 8: Contributions to the nuclear /i — e conversion: (a) long-distance and (b) short-distance. 




Figure 9: LFV branching ratios in the supersymmetric inverse seesaw model of neutrino mass (see text). 

The novel feature present in this model and not in the minimal seesaw is the possibility of enhancing 
Br{fi — > 67) and other tau decays with lepton flavour violation even in the absence of supersymmetry 
in the case where M is low, around TeV or so. In this region of parameters the model also gives rise 
to large estimates for the nuclear — conversion, depicted in Fig. 8. The latter fall within the 
sensitivity of future experiments such as PRISM [85]. Note that LFV happens even in the absence 
of supersymmetry and even in the massless neutrino limit. The allowed rates are hence unsuppressed 
by the smallness of neutrino masses [80, 81, 82, 83, 84]. Finally, for low enough M the corresponding 
quasi-Dirac heavy leptons could be searched directly at accelerators [86, 87]. 

Supersymmetric neutrino mass and collider tests 



17 



10"" 10"^' 10"^^ 

Br(/i ^ ey) 

Figure 10: Correlation between Br{n — > 67) and muon-electron conversion in nuclei from Ref. [79]. 

We now turn to the case of low-scale models of neutrino mass, considered in Sec. 3.2. As an 
example we consider the case of models where supersymmetry is the origin of neutrino mass [53, 58], 
considered in Sec. 3.2.2. A general feature of these models is that, unprotected by any symmetry, 
the lightest supersymmetric particle (LSP) is unstable. In order to reproduce the masses indicated 
by current neutrino oscillation data, typically the LSP is expected to decay inside the detector. More 
strikingly, its decay properties correlate with neutrino mixing angles. For example, if the LSP is the 
lightest neutralino, it is expected to have the same decay rate into muons and taus, since the observed 
atmospheric angle, as we will see below, is close to 7r/4 [88, 89, 90]. This opens the tantalizing possibility 




510^10^ 0.5 1 0.1 0.5 1 5 



Figure 11: LSP decays trace the atmospheric mixing angle [88]. 

of testing neutrino mixing at high energy accelerators, like the "Large Hadron Collider" (LHC) and 
the "International Linear Collider" (ILC) and constitutes a smoking gun signature of this proposal 
that for sure will be tested. This possibility also illustrates the complementarity of accelerator and 
non-accelerator approaches in elementary particle physics. 
Majorons and B-L gauge bosons 

Before closing let us also mention that there could be a dynamical "tracer" of the mass-generation 
mechanism. In the case of the "1-2-3" majoron seesaw considered in Sec. 3.1.1 the existence of the 
Goldstone boson brings new interactions of neutrinos with majoron, but the rates are probably too 
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small to be cosmologically relevant. There can also be a new neutrino decay mode, 1/3 — > + J, 
involving the emission of the majoron, denoted J [8] . However, the strong suppression of the majoron 
emitting decay rates obtained in Ref. [8] and the smallness of neutrino masses that follows from current 
laboratory experiments, as well as cosmology, indicates that such neutrino decays are purely academic. 
However, there is an important exception. In models with family-dependent lepton numbers the decay 
is enhanced even in vacuum [91, 92]. Alternatively, if the neutrino decays in high density media, 
hke supernovae, characterized by huge matter effects, then they could lead to detectable signals in 
underground water Cerenkov experiments [65] even in the simplest majoron schemes. Conversely, in 
the case of seesaw mechanisms with gauged B-L symmetry, if the corresponding B-L scale is low, as 
in the model discussed in Sec. 3.1.4, there will exist a light new neutral gauge boson, Z' that could be 
detected in searches for Drell-Yan processes at the LHC. 

4 The general structure of the lepton mixing matrix 

In any gauge theory in order to identify physical particles one must diagonalize all relevant mass ma- 
trices, which typically result from gauge symmetry breaking. Mechanisms giving mass to neutrinos 
generally imply the need for new interactions whose Yukawa couplings (like Y^^ will coexist with that 
of the charged leptons, Yi. The lepton mixing matrix V follows from a mismatch between the diagonal- 
ization of the charged lepton mass matrix and that of the neutrino mass matrix, which involve these 
Yukawas. This is similar to the way the Cabibbo-Kobayashi-Maskawa matrix arises in the quark sector, 
namely from a mismatch between up and down-type Yukawa couplings. Hence, like quarks, massive 
neutrinos will generally mix. The structure of this mixing is not generally predicted from first principles. 
Whatever the ultimate high energy gauge theory may be it must be broken to the SM at low scales, so 
one should characterize the structure of the lepton mixing matrix in terms of the 5'C/(3)(8)5'C/(2)(g)C/(l) 
structure. 

4.1 Quark- like lepton mixing matrix 

Prom the start, in the absence of a fundamental theory of flavour, the lepton mixing matrix of massive 
Dirac neutrinos is a unitary matrix V which can always be parametrized as 

n 

where 

n 

Wo (7) = expi(^7„^^) 

a=l 

is a diagonal unitary matrix described by n — 1 real parameters 7^. (By choosing an overall relative 
phase between charged leptons and Dirac neutrinos we can take V as unimodular, i. e. det ^ = 1, so 
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that the phases in the "Cartan" matrix loq obey Y12=i 7a — On the other hand 

n 

uJabiVab) = exp J^i'^abAl - rjl,,A^) 

a=l 



is a complex rotation in ab with parameter rjab — \Vab\^^P i^ab- For example, 

/ 

V J 



Cl2 




0... 


\ 




Cl2 


0.. 










1... 





(41) 



However, in order to obtain the parametrization of this matrix, as presented in Ref. [7], one must take 
into account that its structure can be further simplified by taking into account that, once the charged 
leptons and Dirac neutrino mass matrices are diagonal, one can still rcphasc the corresponding fields 
by u!o{a) and uJoil — ct), respectively, keeping invariant the form of the free Lagrangean. This results 
in the form 

n 

V = uJo{a)YluJij{r]ij) ujl{a). (42) 

i<j 

There are n — 1 a-values associated to Dirac neutrino phase redefinitions which we are still free to 
choose. Using the conjugation property 

ujo{a)uJab{\riab\exp iOab) iol{a) = ujab[\riab\exp i{aa + 9ab - ttb)] (43) 

we arrive at the final Dirac lepton mixing matrix which is, of course, identical in form to that describing 
quark mixing. It involves a set of 

n{n — l)/2 mixing angles 9ij and n{n — l)/2 — (n — 1) CP phases . (44) 

One sees that (n — 1) phases were eliminated by rephasing, a possibility that would be absent in 
the case of Majorana neutrinos. This is the parametrization as originally given in [7] , with unspecified 
factor ordering. It needs to be supplemented only by an ordering prescription, for this we refer the 
reader to the PDG choice [93]. 

In summary, if neutrinos masses were added a la Dirac their charged current weak interaction would 
have exactly the same structure as that of quarks. Prom Eq. (44) for the case n — 3 there are 3 angles 
and precisely one leptonic CP violating phase, just as in the Kobayashi-Maskawa matrix describing 
quark mixing. 

However the imposition of lepton number conservation is ad hoc in a gauge theory and hence 
neutrinos are expected to be Majorana ^. Note that the argument in favour of neutrinos being Majorana 
has nothing to do with the neutrino mass generation mechanism, it is much more basic. 



^The same happens for electrically neutral supersymmetric fermions such as the gravitino, gluino and neutralinos. 
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4.2 Majorana case: unitary approximation 



Here we consider the form of the lepton mixing matrix in models where neutrino masses arise in the 
absence of right-handed neutrinos. These include, for example, the models in Sec. 3.2.1 and constitute 
a good approximation for high-scale seesaw models in Sec. 3.1. 

For n generations of Majorana neutrinos the lepton mixing matrix has exactly the same form given 
in Eq. (42). The difference insofar as the structure of lepton mixing is concerned come about since in the 
case of Majorana neutrinos the mass terms are manifestly not invariant under rephasings of the neutrino 
fields. As a result, the parameters a in Eq. (42) can not be used to eliminate n — 1 Majorana phases 
as we just did in Sec. 4.1. Consequently these are additional sources of CP violation in the currents of 
gauge theories with Majorana neutrinos. They are sometimes called "Majorana phases". They already 
exist in a theory with just two generations of Majorana neutrinos, n = 2, which is described by 




Cl2 



(45) 



where 0i2 is the Majorana phase. Such "Majorana" -type CP phase is, in a sense, mathematically more 
"fundamental" than the Dirac phase whose existence, as we just saw, requires three generations at least. 
For the case of three neutrinos the lepton mixing matrix can be parametrized as (n = 3) [7] 



(46) 



where each factor in the product of the a;'s is effectively 2x2, characterized by an angle and a CP 
phase. Two of the three angles are involved in solar and atmospheric oscillations, so we set 9i2 = ^sol 
and 023 = ^ATM- The last angle in the three-neutrino leptonic mixing matrix is ^13, 



V 



Cl3 

1 
-^"^"sia 





Cl3 



(47) 



Such symmetrical parametrization of the lepton mixing matrix, K can be written as: 



K ^ 



C12C13 



Sl2Cl3e 



l(pl2 



l<Pl3 



-Sl2C23e-''^i2 - Ci2Si3S23e*('^^^-'^i3) C12C23 - Si2Sl3S23e'^'l'^^+'^^^-'f>'-^^ CiiS23e''^^' 



512^236 



-i{4>12+4'23) 



Cl2Sl3C23e 



-Cl2S23e 



-i<t>23 



Sl2Sl3C23e 



i(</>12— 013) 



C13C23 



(48) 



All three CP violating phases are physical [25]: <^i2,</'23 and (pi^. Even though the parametrization 
is fully "symmetric" there is a basic difference between Dirac and Majorana phases. The "invariant" 
combination 6 = 0i2 + 023 — 0i3 corresponding to the "Dirac phase" exists only beyond n = 3. At 
n = 3 a single phase (say ^is) may be taken to be non-zero. This is the phase that corresponds to 
the Dirac phase present in the quark sector, and affects neutrino oscillations involving three neutrinos. 
The other two phases are associated to the Majorana nature of neutrinos. These already exist for 
n = 2 but show up only in lepton-number violating processes, like neutrinoless double beta decay. They 
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do not affect conventional neutrino oscillations [25, 26, 27] but do enter in the "neutrino oscillation 
thought-experiment" described in [25]. Needless to say the latter process is of conceptual interest, only. 

An important subtlety arises regarding the conditions for CP conservation in gauge theories of 
massive Majorana neutrinos. Unlike the case of Dirac fermions, where CP invariance implies that the 
mixing matrix should be real, in the Majorana case the condition one needs is 

K* = Kri 

where rj — diag(-|-, +,...,,,..) is the signature matrix describing the relative signs of the neutrino mass 
eigenvalues that follow from diagonalizing the relevant Majorana mass matrix, if one insists in using 
real diagonalizing matrices, like Wolfenstein [32]. Consequently the value 0i2 — 7r/2 and 0i2 = are 
both CP conserving. These important signs determine the CP properties of the neutrinos and play a 
crucial role in /3/3oi/ ■ 

Before concluding let us mention that the above parametrization of the lepton mixing matrix was 
originally given in [7], but with unspecified factor ordering. In what follows we tacitly employ the 
ordering prescription now adopted by the PDG [93]. Further discussion on the advantages of this 
"symmetrical parametrization" and additional references can be found in Ref . [94] . 

4.3 General seesaw- type lepton mixing matrix 

What we now present is a brief summary of the original systematic study of the effective form of the 
lepton mixing matrix originally given in [7], applicable to any scheme of Majorana neutrino masses 
where isosinglet and isodoublet mass terms coexist. In addition to the existence of Majorana phases, 
which do not require the existence of SU{3) (8) SU{2) (g) U{1) singlets, one has doublet-singlet mixing 
parameters, in general complex. As a result one finds that leptonic mixing as well as CP violation may 
take place even in the massless neutrino hmit [81, 82]. 

The most general effective model is described by (n, m), n being the number of SU{S)®SU{2)®U{1) 
isodoublet and m the number of SU{3) SU{2) (8) C/(l) isosinglet leptons. In this Section we assume 
m 7^ 0, the (n, 0) has just been considered in Sec. 4.2. The two-component SU{3) ® SU{2) ® U{1) 
singlet leptons present in the theory have in general a gauge and Lorentz invariant Majorana mass 
term, breaking total lepton number symmetry. Their number, m, is completely arbitrary, as SU (3) ® 
SU{2) ® ?7(1) singlets carry no anomaly. 

The resulting structure of the weak currents can be substantially more complex, since the heavy 
isosinglets will now mix with the ordinary SU{2) doublet neutrinos in the charged current weak inter- 
action. As a result, the mixing matrix describing the charged leptonic weak interaction is a rectangular 
matrix, called K. 

Typically the charged weak interactions of the light (mass-eigcnstate) neutrinos in {n, m) models 
arc effectively described by a mixing matrix which is non-unitary. For example, one sees that the 
coupling of a given light neutrino to the corresponding charged lepton is decreased by a certain factor. 
The existence of these neutral heavy leptons could therefore be inferred from low energy weak decay 
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processes, where the neutrinos that can be kinematically produced are only the hght ones. There are 
constraints on the strength of the such mixing matrix elements that follow from low energy weak decay 
measurements as well as from LEP. 

An explicit parametrization of the weak charged current mixing matrix K that covers the most 
general situation present in these (n, m) models has also been given in Ref. [7] so that here we only 
highlight the main points. It involves in general 

n(n + 2m- l)/2 (49) 

mixing angles 9ij and 

n(n + 2m-l)/2 (50) 

CP violating phases 

This number far exceeds the corresponding number of parameters describing the charged current 
weak interaction of quarks Eq. (42). The reason is twofold: (i) neutrinos are Majorana particles, their 
mass terms are not invariant under rephasings, and (ii) the isodoublet neutrinos in general mix with 
the SU{?>) ® SU{2) (g) C/(l) singlets and, in addition, so CP may also be violated in this mixing. As a 
result, there are far less CP phases that can be eliminated by field redefinitions. They may play a role 
in leptogenesis as well as neutrino oscillations. 

Another important feature which arises in any theory based on SU{3) SU{2) <Si U{1) where 
isosinglet and isodoublet lepton mass terms coexist is that the leptonic neutral current is non-trivial [7] : 
there are non diagonal couplings of the Z to the mass-eigenstate neutrinos. They are expressed as a 
projective Hermitian matrix 

P^K^K 

This contrasts with the neutral current couplings of mass-eigcnstatc neutrinos in theories where there 
are no isosinglet neutrinos , i.e., m = 0. In that case, just as in the case of Dirac neutrinos, the neutral 
current couplings of mass-eigcnstatc neutrinos is diagonal. 

Before closing this section we mention explicitly that the (3, 3) seesaw model is characterized, from 
Eqs. (49) and (50) by 12 mixing angles and 12 CP phases (both Dirac and Majorana-type) [7]. 

Before we close, note that, in a scheme with m < n, n — m neutrinos will remain massless, while 
2m neutrinos will acquire Majorana masses, m light and m heavy. For example, in a model with n — 3 
and m = 1 one has one light and one heavy Majorana neutrino, in addition to the two massless ones. 
In this case clearly there will be less parameters than present in a model with m = n. Note also that 
for m > n, the case m — 2n corresponds to the models considered in Sees. 3.1.4 and 3.2.3. 

4.4 Leptonic CP violation and leptogenesis 

An interesting cosmological implication of leptonic CP violation is that it opens an attractive possibility 
of accounting for the matter-antimatter asymmetry in the Universe within the so-called thermal lepto- 
genesis mechanism [95, 96]. In this picture the heavy "right-handed" neutrinos of the seesaw mechanism 
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play a crucial role. They decay out of equilibrium generating a lepton asymmetry, through diagrams 
in Fig. 12. This lepton (or B-L) asymmetry gets converted, through sphaleron processes, into a baryon 



Figure 12: Diagrams contributing to heavy neutrino decays whose interference leads to leptogenesis. 

asymmetry [97]. One of the crucial ingredients is CP violation in the lepton sector. 

In the framework of a supersymmetric seesaw scheme the high temperature needed for leptogenesis 
leads to an overproduction of gravitinos, which destroys the standard predictions of Big Bang Nucle- 



stable, decaying during or after BBN. Their rate of production can be so large that subsequent gravitino 
decays completely change the standard BBN scenario. To prevent such "gravitino crisis" one requires 
an upper bound on the reheating temperature Tr after inflation, since the abundance of gravitinos is 
proportional to the reheating temperature. A recent detailed analysis derived a stringent upper bound 
Tr ^ 10^ GeV when the gravitino decay has hadronic modes [98]. This upper bound is in conflict 
with the temperature required for leptogenesis, Tji> 2x 10^ GeV [99]. Therefore, thermal leptogenesis 
in the usual hierarchical case seems difficult to reconcile with low energy supersymmetry if gravitino 
masses lie in the range suggested by the simplest minimal supergravity models. Their required mass is 
typically too large in order for them to be produced after inflation, implying that the minimal type I 
supersymmetric seesaw schemes may be in trouble ^. 

One way to cure this inconsistency has been suggested in Ref. [101] is to add a singlet E to the 
supersymmetric SO{10) seesaw model discussed in [48]. Leptogenesis can occur at the TeV scale through 
the decay of E, thereby avoiding the gravitino crisis. Washout of the asymmetry is suppressed by the 
absence of direct couphngs of E to leptons. The leptogenesis parameter region is illustrated in Fig. 13 
and the reader is addressed to the original paper for details. Here we simply comment that in this model 
successful leptogenesis can occur for Ms = 1 TeV and low vr — 10 TeV. As an example we consider 
the E decay asymmetries ei that can be generated in this way. The contours of the maximum ei values 
obtained in this model using the Fritzsch texture for the quark masses are illustrated in Fig. 14. The 
left and right panels correspond to two alternative choices of the relevant model parameters, explained 
in Rcf. [102]. Taking into account washout effects, acceptable values of the baryon asymmetry O(10^^°) 
require ei above (9(10~^) or so. From Fig. 14 one sees that this is clearly possible to achieve within this 
model even if the only source of CP violation is that provided by the three-flavour Dirac phase that 
enters in neutrino oscillations. 

loophole is the possibility of resonant leptogenesis [100]. 




osynthesis (BBN). In minimal supergravity models, with 7713/2 ~ 100 GeV to 10 TeV gravitinos are not 
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10° 10^ 10-^ lO'' 10^ VR[TeV] 

Figure 13: Parameter space in low-scale leptogenesis model of Ref. [101]. 




An alternative viable leptogenesis mechanism based on a modified type-I seesaw was suggested in 
Ref. [103]. It consists in adding a small R-parity violating XiV'^iH^Hct term in the superpotential, where 
v'^i are right-handed neutrino supermultiplets. One can show that in the presence of this term, the 
produced lepton-antilepton asymmetry can be enhanced. 

Last, but not least, leptogenesis may also be induced by the type-II seesaw mechanism, with or 
without supersymmetry. In the non-supersymmetric case this requires the addition of two heavy Higgs 
scalar triplets [104], while in supersymmetry even with a minimal triplet content, leptogenesis can be 
naturally accommodated thanks to the resonant interference between superpotential and soft super- 
symmetry breaking terms [105]. 
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5 Neutrino oscillations 



5.1 The formalism 

We start with a very general discussion aimed at fixing the notation used in this review. As we have 
already seen in Sec. 4.3, in general models of neutrino mass the lepton mixing matrix K contain both 
Dirac and Majorana-type CP phases and is in general rectangular as it couples the charged leptons 
also to the heavy (mainly isosinglet) neutrinos postulated, c. g. in typc-1 seesaw models, in order to 
produce neutrino masses [7] . Such states are too heavy to participate in neutrino oscillations which are 
effectively described by a non-unitary mixing matrix. Such deviations from unitarity arc the origin of 
gauge-induced neutrino non-standard interactions and may, in some cases, be sizeable, sec Sec. 9. In 
the discussion of neutrino oscillations that we give in Sees. 5, 6 and 7 we will tacitly assume that K is 
strictly unitary, so that the three active neutrinos are mixed as follows, 



(51) 



where we have now denoted K by U, to highlight that U^U — UW — 1. Here u^l (« = 
describe the left handed neutrino fields with definite flavor whereas {i — 1,2,3) describe the fields 
with definite masses. Here, the matrix U is the leptonic analogue of the quark mixing matrix [1]. In 
this review, whenever it is necessary to use the explicit parametrization, we use the following standard 
parametrization [93], 
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(52) 



(53) 



where = cos^y, = sin 9ij and S is the CP violating phase. This form is the same as Eq. (46) [7], 
taking 0i2 = = 023 in the "invariant" combination 5 = 0i2 + 023 — 0i3 that corresponds to the "Dirac 
phase" relevant for neutrino oscillations. 

Eq. (51) implies that in terms of state vectors {a — e,iJ,,T) for flavor and = 1,2,3) for 

mass eigenstates, they are related by U* instead of U as (see e.g. [106]), 



(54) 



i=l 



where for simplicity, we omitted the indices L which indicates the left handed chirality. Then it is 
straightforward to compute the Ua I'p oscillation probability, which is given, for ultra-relativistic 
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neutrinos, by, 



i>j 



2 I 



i>j ^ / 



(55) 



where E is the neutrino energy, L is the distance traveled by neutrino, and Am^j = mf — {rrii being 
mass eigenvalues) are the mass squared differences. Here 3? and 9 denote real and imaginary parts. 

Let us now describe the neutrino evolution equation in matter. It is given in terms of flavor eigen- 
states as. 



. d 

dx 



where (q^ = A*, t) is the amplitude for the a-fiavor. The Hamiltonian matrix H is given by 
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(57) 



where x is the position along the neutrino trajectory and V{x) is the matter potential which is given, 
for unpolarized medium ^, as 

V(x) = V2GFNe{x) , (58) 

where is the Fermi constant, Ne{x) is the electron number density at x. The Hamiltonian matrix 
in Eq. (57) can be replace by. 



H{x) = C/diag 



0, 
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2E ' 2E 



[/^ + diag[y(x),0,0]. 



(59) 



by re- phasing all of the neutrino flavors by exp[— imix/(2£')]. For anti-neutrinos, the same Eq. holds 
with the change V{x) — > —V{x) and U ^ U* (or equivalently 5 — > —5). 

For a given x, Hamiltonian can be diagonalized by using the effective mixing matrix in matter U (N), 



H{x) = C/(A^)diag 



0, 



Amii(iV) Aml^{N) 



U\N), 



2E ' 2E 

where AmfAN) are the effective mass squared differences in matter (see Eq. 88). 



(60) 



^Neutrino evolution in polarized media was treated in [107]. 
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5.2 The current data 



Current neutrino oscillation data have no sensitivity to CP violation. Thus we neglect all phases in the 
analysis and take, moreover, the simplest unitary 3-by-3 form of the lepton mixing matrix in Eq. (46). 
In this approximation oscillations depend on the three mixing parameters sin^ ^^12, sin^ ^^23i sin^ ^^i3 and 
on the two mass-squared splittings Am'^Q^ = Amli = m\ — m\ and Am^^^j = A.m\^ = m\ — m\ 
characterizing solar and atmospheric neutrinos. The hierarchy Am^Q^ ^ ^"^atm implies that one can 
set Am^QL = 0, to a good approximation, in the analysis of atmospheric and accelerator data. Similarly, 
one can set Am^^^, to infinity in the analysis of solar and KamLAND data. Apart from the data already 
mentioned, the analysis also includes the constraints from "negative" searches at reactor experiments. 

5.2.1 Solair and reactor data 

The solar neutrino data includes the rates of the chlorine experiment (2.56 ± 0.16 ± 0.16 SNU), the 
gallium results of SAGE (66.9 3;^ SNU) and GALLEX/GNO (69.3 ± 4.1 ± 3.6 SNU), as well as 
the 1496-day Super-K data (44 bins: 8 energy bins, 6 of which are further divided into 7 zenith angle 
bins). The SNO sample includes the 2002 spectral day/night data (17 energy bins for each day and 
night period) [3] and the most recent data (391-day data) from the salt phase in the form of the neutral 
current (NC), charged current (CC) and elastic scattering (ES) fluxes [108]. The analysis includes 
both statistical errors, as well as systematic uncertainties such as those of the eight solar neutrino 
fluxes. Taking into account new radiative opacities, Bahcall et al (for references see Appendix C in hep- 
ph/0405172-v5) obtain new solar neutrino fluxes, neutrino production distributions and solar density 
proflle, included in the analysis. 

Reactor anti-neutrinos from a network of power stations in Japan are detected by the KamLAND 
collaboration at the Kamiokande site through the process Pe+P ~^ e'^-\-n, where the delayed coincidence 
of the prompt energy from the positron and a characteristic gamma from the neutron capture allows 
an efficient reduction of backgrounds. Most of the incident V^s come from nuclear plants at distances 
of 80 — 350 km from the detector, far enough to probe large mixing angle (LMA) oscillations. To avoid 
large uncertainties associated with geo-neutrinos an energy cut at 2.6 MeV prompt energy is applied 
for the oscillation analysis. 

The flrst KamLAND data corresponding to a 162 ton-year exposure gave 54 anti-neutrino events in 
the flnal sample, after cuts, whereas 86.8 ± 5.6 events are predicted for no oscillations with 0.95 ± 0.99 
background events [109]. This is consistent with the no-disappearance hypothesis at less than 0.05% 
probability, giving the first terrestrial confirmation of oscillations with ArngQ^- Additional KamLAND 
data with a larger fiducial volume of the detector corresponding to an 766.3 ton-year exposure have 
been presented in [4]. In total 258 events have been observed, versus 356.2±23.7 reactor neutrino events 
expected in the case of no disappearance and 7.5 ± 1.3 background events. This leads to a confidence 
level of 99.995% for disappearance, in addition to evidence for spectral distortion consistent with 
oscillations. 



28 



Finally, recent data from the KamLAND experiment have been presented at the 10th International 
Conference on Topics in Astroparticle and Underground Physics, TAUP 2007 [110]. These data corre- 
spond to a total exposure of 2881 ton-year. They provide a better measurement of the solar neutrino 
oscillation parameters, especially the "solar" mass-squared splitting. This is due to the reduction of 
systematic uncertainties thanks to the full volume calibration. These data correspond to a total expo- 
sure of 2881 ton-year, almost 4 times larger than 2004 data. They provide a very precise measurement 
of the solar neutrino oscillation parameters, mainly the mass splitting. This is due to the reduction of 
systematic undertainties thanks to the full volume calibration. 

Various systematic errors associated to the neutrino fluxes, backgrounds, reactor fuel composition 
and individual reactor powers, small matter effects, and improved Pg flux parametrization are included 
in the analysis [12]. Assuming CPT invariance one can directly compare the information obtained from 
solar neutrino experiments with the KamLAND reactor results. 

One finds that a strong evidence for spectral distortion in the KamLAND data, leading to a much 
improved AmgoL determination, substantially reducing the allowed region of oscillation parameters. Al- 
together the KamLAND data single out the LMA solution from the previous "zoo" of alternatives [111]. 
As discussed in Sec. 9, more than just cornering the oscillation parameters [112], KamLAND has elimi- 
nated all previously viable non-oscillation solutions [113] playing a key role in the resolution of the solar 
neutrino problem. 

Last, but not least, the Borexino collaboration has also recently presented their first data [114]. 
While they provide the first real time detection of ^Be solar neutrinos, and an important confirmation 
of the Standard solar model and the large mixing oscillations, these data currently do not affect the 
determination of neutrino oscillation parameters. 

5.2.2 Atmospheric and accelerator data 

The first evidence for neutrino conversions was the zenith angle dependence of the /x-like atmospheric 
neutrino data from the Super-K experiment in 1998, an effect also seen in other atmospheric neutrino 
experiments. At that time there were equally good non-oscillation solutions, involving non-standard 
neutrino interactions [115]. Thanks to the accumulation of up-going muon data, and the observation of 
the dip in the L/E distribution of the atmospheric survival probability, the signature for atmospheric 
neutrino oscillations has now become convincing. We have used the Super-K charged-current atmo- 
spheric neutrino events, with the e-like and /x-like data samples of sub- and multi-GeV contained events 
grouped into 10 zenith- angle bins, with 5 angular bins of stopping muons and 10 through-going bins of 
up-going muons [5]. Multi-ring // and neutral-current events and Vj- appearance are not used, since an 
efficient Monte-Carlo simulation of these data would require further details of the Super-K experiment, 
in particular of the way the neutral-current signal is extracted from the data. As far as atmospheric neu- 
trino fluxes are concerned, the analysis of [12] employs state-of-the-art three-dimensional calculations 
given in [17]. 
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The disappearance of i/^'s over a long-baseline probing the same Am^ region relevant for atmo- 
spheric neutrinos is now available from accelerator neutrino oscillation experiments, such as the KEK 
to Kamioka (K2K) neutrino oscillation experiment and the MINOS Experiment using the NuMI Beam- 
line facility at Fermilab. 

Neutrinos produced by a 12 GeV proton beam from the KEK proton synchrotron consist of 98% 
muon neutrinos with a mean energy of 1.3 GeV. The beam is controlled by a near detector 300 m away 
from the proton target. Comparing these near detector data with the content of the beam observed 
by the Super-K detector at a distance of 250 km gives information on neutrino oscillations. 

The K2K collaboration has pubhshed details of the analysis of their full data sample (K2K-I and 
K2K-II) [119]. The data have been taken in the period from June 1999 to November 2004 and correspond 
to 0.922 X 10^° p.o.t. Without oscillations 158l'l8;6 events are expected whereas only 112 events have 
been observed. Out of these, 58 events are single-ring events where the reconstruction of the neutrino 
energy is possible. These events have been used in order to perform a spectral analysis of the K2K 
data, as described in Ref. [12]. 

The first MINOS results were released in 2006. MINOS is a long-baseline experiment that searches 
for disappearance in a neutrino beam with a mean energy of 3 GeV produced at Fermilab. It consists 
of a near detector, located at 1 km from the neutrino source and a far detector located at the Soudan 
Mine, at 735 km from Fermilab. First data corresponding to 0.93 x 10^° p.o.t. [117], have comparable 
weight as the final K2K data sample. In the absence of oscillations 177 ± 11 events with £' < 10 GcV 
are expected, whereas 92 have been observed, which provides a 5.0(T evidence for disappearance. New 
experimental data have been released by the MINOS Collaboration. These have been collected from 
June 2006 to July 2007 (Run-IIa), and they have been analyzed together with the first data sample 
(Run-I), with a total exposure of 2.5x10^° p.o.t. In total, 563 events have been observed at the far 
detector, while 738±30 events were expected for no oscillation. 

One finds that the values of the oscillation parameters from the z/^ disappearance results at MINOS 
arc consistent with the ones from K2K, as well as from Supcr-K atmospheric data, providing strong 
terrestrial confirmation of oscillations with Am^^j^j with accelerator neutrinos. Statistics of the current 
data sample still does not strongly constrain the mixing angle. However, although the determination 
of sin ^ATM is completely dominated by atmospheric data, K2K and MINOS data give an important 
restriction on the allowed Am^,j,^ values [12]. They constrain Am^^^^, from below, which is important for 
further future long-baseline experiments, since their sensitivities arc drastically affected if Am^^^j lies 
in the lower part of the 3a range indicated by current atmospheric data. Note that there is currently 
no sensitivity to the sign of the atmospheric mass-squared splitting, this being one of the challenges for 
upcoming experiments. 

5.3 Status of three-neutrino oscillations 

The three-neutrino oscillation parameters that follow from the global analysis of neutrino oscillation 
data have been reviewed by Maltoni et al in Ref. [12]. An update of the results of this analysis 
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Figure 15: Current regions of neutrino oscillation parameters allowed by the world's neutrino oscillation 
data at 90%, 95%, 99%, and C.L. for 2 d.o.f. from Ref. [12]. Shaded regions include 2007 data, 
while the empty regions correspond to results before the latest update. 



can be found in the arXiv as version 6 of hep-ph/0405172. This update includes the latest relevant 
experimental data, as well as the latest update of the SSM of Bahcall et al (for references see Appendix 
C in hep-ph/0405172-v6). The results are summarized in Figs. 15 and 16 as well as in the Table. 

The plots shown in Fig. 15 are two-dimensional projections of the allowed regions in the five- 
dimensional parameter space. The left panel in Fig. 15 gives the allowed regions of the "solar" neutrino 
oscillation parameters (sin^6'i2, Amgi). One sees the important role of the 2007 KamLAND data in 
providing an improved determination of both parameters, especially the mass squared splitting, now 
displayed on a linear scale. In the right panel of Fig. 15 we give the allowed regions of the "atmospheric" 
parameters (sin^ 6'23, Am|^) also on a linear scale. In addition to a confirmation of oscillations with 
Am^^^ with accelerator neutrinos, these data provide a better determination of Am^^^,, with little 
impact on the other parameters. Indeed, one can see from the figure that the inclusion of the MINOS 
data is crucial in the determination of Am^^^j. We see how accelerator data are starting to play the 
main role in the determination of the "atmospheric" splitting, a trend which will become stronger in 
the future. The best fit values and the allowed 3cr ranges of the oscillation parameters from the global 
data are summarized in the Table. 

The left panel in Fig. 16 gives the parameter a, namely the ratio of solar over atmospheric splittings, 
as determined from the global analysis. The right panel in Fig. 16 illustrates how the bound on 
sin^ ^13 emerges from the interplay of different data samples. The plot shows the upper bound on sin^ ^13 
as a function of /S.m\^^^ from CHOOZ data alone and compares to the bound obtained from an analysis 
including also solar and reactor neutrino data. One sees that, although for larger Am^^^^ values the 
bound on sin^ 6'i3 is dominated by CHOOZ, this bound deteriorates quickly as Am^^^j decreases, so that 
for Am^^^j ^ 2 X 10~'^eV^ the solar and KamLAND data become relevant. 
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Figure 16: Current determination of a = Am^Qj^/Am^^^j and bound on sin^ Oi^ from the world's neutrino 
oscillation data, from Ref. [12]. 



parameter 


best fit 


2a 


3a 


Amj^ [lO-^eV^] 


7.6 


7.3-8.1 


7.1-8.3 


Amii [lO-^eV^] 


2.4 


2.1-2.7 


2.0-2.8 


sin^ 6*12 


0.32 


0.28-0.37 


0.26-0.40 


sin^ 023 


0.50 


0.38-0.63 


0.34-0.67 




0.007 


< 0.033 


< 0.050 



Table I: Best-fit values, 2a and 3a intervals (1 d.o.f.) for the three-flavour neutrino oscillation 
parameters from global data including solar, atmospheric, reactor (KamLAND and CHOOZ) and 
accelerator (K2K and MINOS) experiments. 

In summary, we find the following bounds at 90% C.L. {3a) for 1 d.o.f.: 

{ 0.051 (0.084) (solar+KamLAND) 
0.028 (0.059) (CHOOZ+atmospheric+K2K+MINOS) (61) 
0.028 (0.050) (global data) 

5.4 Predicting neutrino oscillation parameters 

As we saw in Sec. 5 only five of the basic parameters of the lepton sector are currently probed in 
neutrino oscihation studies [2, 3, 4, 5, 6] [14, 15]: the angle 6*12 and the splitting Am^Q^^, which are 
determined from solar and KamLAND data, together with the angle 623 and the corresponding mass 
squared splitting Am^^,^ determined by atmospheric and accelerator data. 
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Out of the five neutrino oscillation parameters in the three-neutrino leptonic mixing matrix, two 
play a key role for future CP violation studies in neutrino oscillation experiments: the third "small" 
angle ^13, together with the ratio of solar to atmospheric splittings, derived from the oscillation analysis. 
These are given in Fig. 16. 

As we have seen, current neutrino data point towards a well defined pattern of neutrino mixing 
angles, quite distinct from that of quarks. The data seem to indicate an intriguing complementarity 
between the angles that characterize the quark and lepton mixing matrices [120, 121, 122, 123]. It is 
not clear whether this numerological coincidence has a deeper meaning. 

Some "flavour-blind" gauge models have the interesting feature of "accidentally" predicting a zero 
neutrino mass. For example, due to the anti-symmetry of the one of the Yukawa coupling matrices, the 
radiative models [50, 51] considered in Sec. 3.2.1 predict that one of the neutrinos is massless. Similarly 
for the model considered in Ref. [30]. These are rather exceptional examples and the predictions 
obtained are incomplete and/or problematic, such as those of the simplest Zee model [50] which also 
include a nearly maximal solar mixing angle, now in conflict with observation. This illustrates that, 
typically, gauge symmetry by itself is not sufficient to predict masses and mixings, neither for quarks, 
nor for leptons. 
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Figure 17: Near maximal CP violation in neutrino oscillations, as predicted in Ref. [133]. 

There has been a rush of papers attempting to understand the values of the leptonic mixing angles 
from underlying symmetries at a fundamental level. For example one can account for the maximum 
mixing in atmospheric neutrino oscillation, as well as vanishing of the 613 

023 = 7r/4 and ^13 = (62) 

in terms of a mu-tau symmetry, under which the neutrino mass matrix is invariant under the inter- 
change of second and third generation neutrinos [124]. Other models with Z2 symmetries have been 
considered [125]. For example, in [126] one such symmetry has been considered leading to a maximal at- 
mospheric mixing angle 623, its breaking inducing a nonzero value of ^13 in a way that strongly depends 
on the neutrino mass hierarchy. 
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Figure 18: Magnitude of CP violation invariant in neutrino oscillations, as predicted in Ref. [133]. 
More generally, it has been argued that the form of the mixing at high energies could be [127] 
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Such Harrison-Perkins-Scott pattern of lepton mixing [128] could result from some kind of flavour 
symmetry, valid at a very high energy scale where the dimension-five neutrino mass operator arises. 

One idea is that neutrino masses arise from a common seed at some "neutrino mass unification" 
scale Mx [129], very similar the merging of the gauge coupling constants at high energies due to 
supersymmetry [130]. Unfortunately in its simplest (CP conserving) form this very simple theoretical 
ansatz is inconsistent with the significantly non-maximal value of the solar mixing angle 612 inferred 
from current data [12]. 

A more satisfactory and fully viable alternative realization of the "neutrino mass unification" idea 
employs an A4 fiavour symmetry in the context of a seesaw scheme [70]. Starting from three- fold 
degeneracy of the neutrino masses at a high energy scale, a viable low energy neutrino mass matrix can 
indeed be obtained. The model predicts maximal atmospheric angle and vanishing 613, as in Eq. (62). 
Moreover, if CP is violated ^13 becomes arbitrary but the Dirac CP violation phase is maximal [131]. 
The solar angle 6^12 is unpredicted. However, one expects it to be large, 

9i2 = O (1). 

(There have been variant realizations of the symmetry that enable one to predict the solar angle, 
see, for example Ref. [132]). Within such ^4 fiavour symmetry seesaw scheme one can show that the 
lepton and slepton mixings are intimately related. It was shown that the resulting slepton spectrum 
must necessarily include at least one mass eigenstate below 200 GeV, which can be produced at the 
LHC. The prediction for the absolute Majorana neutrino mass scale 

mo > 0.3eV 
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ensures that the model will be tested by future cosmological tests and (3(3qi, searches. Rates for lepton 
flavour violating processes ij — > £j + 7 in the range of sensitivity of current experiments are typical in 
the model, with BR(// 67) ^ 10"^^ and the lower bound BR(t /X7) > 10~^. 

A variant A4 flavour symmetry seesaw model has been suggested [133] which predicts a direct 
correlation between the expected magnitude of CP violation in neutrino oscillations and the value of 
sin^ ^13. The model leads to nearly maximal leptonic CP violation in neutrino oscillations through the 
CP phase S. These predictions are illustrated in Figs. 17 and 18. For a discussion of various schemes 
to derive the neutrino mixing matrix using the tetrahedral group A4 was given in Ref . [134] . 

Also discrete flavour groups of the series have been used. For example, seesaw models for the 
lepton sector using and 5*3 family symmetries have been suggested in [135]. The model predicts 
a normal hierarchy neutrino mass spectrum with the mixing angle ^^13 = and an unpredicted solar 
mixing angle 612- It employs an enlarged Higgs scalar sector, so that the atmospheric mixing angle 623 
is given as a ratio of Higgs scalar vevs, and it is maximal if the full Lagrangian is Z}4-invariant. The 
deviation of 623 from 7r/4 is governed by the strength of the soft breaking of the D4 symmetry. 

Other non-Abclian finite subgroups of SU(3) have also been used as family symmetries and shown 
to generate tri-bimaximal mixing in the neutrino sector, while allowing for quark and charged lepton 
hierarchies [136]. Many attempts at predicting lepton mixing angles employ the Harrison-Perkins-Scott 
mixing pattern in Eq. (63) at some high energy scale [137]. Typically one must correct predictions made 
at a high scale by renormalization group evolution [138]. In this connection the idea of tri-bimaximal 
neutrino mixing has also been considered in the context of discrete symmetries in models with extra 
dimensions [139]. 

In short, progress made in the experimental determination of neutrino parameters has already ruled 
out many of the proposed models. For example, many interesting attempts have failed simply because 
the solar mixing angle has been shown to be large but significantly non-maximal [12]. 

Note that understanding the observed pattern of masses and mixings in the quark and lepton 
sector separately in the context of a non-unified model by appealing to suitable fiavour symmetries is 
a relatively feasible challenge. However, in a unified model where quarks and leptons sit in the same 
multiplets, there is a tendency to correlate the corresponding mixing angles. In a "generic" fiavour-blind 
unified model the lepton and quark mixing angles can always be reconciled thanks to the large number 
of parameters involved. For example, we have seen how in general seesaw-type models the leptonic 
mass matrices involves couplings which are absent for the quarks. Even in the simplest type-I seesaw 
mechanism discussed in Sec. 3.1 the B-L violating terms bring in their own independent flavor structure, 
distinct from that of the Dirac mass term which is restricted by the quark sector. Acceptable mixing 
patterns can certainly be accommodated, but not predicted. To achieve some degree of predictivity one 
must appeal to additional symmetries, beyond the gauge symmetry, realized at the unified level and 
suitably chosen so as to restrict the flavour "textures" . This constitutes a rather demanding challenge, 
so far not achieved with full success. Despite some progress on how to obtain successful unifled models 
of flavour [140] using discrete non-abelian groups, we currently lack a fully satisfactory unifled theory 
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of flavour, and it is likely that the "flavour problem" will indeed remain with us for a while. 

6 CP violation in neutrino oscillations 

In gauge theories of neutrino masses, the lepton mixing matrix typically contains a number of CP 
violating phases that may affect neutrino oscillations [7] . In this section, we discuss theoretical aspects 
of CP violation in neutrino oscillations and the possibility of probing it experimentally at future searches. 

It has been recognized for a long time that the Dirac CP violating phase present in the simplest 
thrcc-ncutrino model could in principle be observed in neutrino oscillation experiments [20]. These 
arc likely to be the most promising way to probe directly the Dirac CP phase present in the neutrino 
mixing matrix, unless Oi^ is too small. If CPT is conserved violation of CP implies that of T. In 
this review, we do not consider the possible violation of CPT in the neutrino sector Earlier works 
around 1980 on CP and/or T violation can be found in Refs. [7], [24, 25, 26] and [143, 144]. It has 
now been understood that the difference between oscillation probabilities for neutrino and anti-neutrino 
is proportional to the leptonic analogue of the CP-invariant factor of the quark sector [145]. In the 
last decade, CP violation in neutrino oscillation has received enormous amount of attention in the 
community [146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156, 157]. Here we consider three active 
neutrinos without sterile neutrino [158, 159, 160] which have been invoked in connection with LSND 
and MiniBoone data [10, 11]. 

In this review, in order keep the discussion as much general as possible, we do not show any kind of 
sensitivity plots for speciflc planned experiment or project, but try to keep the discussion at the level 
of neutrino oscillation probabilities. 

6 . 1 Preliminaries 

The simplest measure of CP violation, which is equivalent to T violation if CPT is conserved, would 
be the difference of oscillation probabilities between neutrinos and anti-neutrinos, Pl^a ^z?) ^-nd 
P{'^a ^p), which is given by [143, 144], 

^PuValB = P{Va l^p) - P{i>a ^ = "16 Ja^ siu A12 siu A23 siu A31, (64) 

where we used the notation, Ajj = AmfjL/AE and 

Jaf3 = '^{UalU*2U^iUi32) = ±J, J = 512012523023-513013 sin 5 (65) 

with positive (negative) sign for (anti-)cyclic permutation of the flavor indices e, /i and r. The parameter 
J is the leptonic analogue of the CP-invariant factor for quarks, the unique and phase-convention- 
independent measure for CP violation [145]. 

^See e.g. [141, 142] for such a possibility. 
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Figure 19: The range of probability of finding the a-fiavor in the i-th mass eigenstate as indicated for 
the two different mass hierarchies for the best fit values of the solar and atmospheric mixing parameters, 
sin^ 2^13 = 0.1 as the CP-violating phase, 6 is varied, extracted from [161]. 

Since we do not yet know the value of S and ^13, the current neutrino data give only an upper bound 
on this quantity, as J ^ 0.04. Although we do not yet know whether the mass hierarchy is normal or 
inverted this will not cause an ambiguity in the vacuum AP, as it does not depend on the hierarchy ^. 

Where we can expect to observe the effect of CP violation? In contrast with CP violation induced 
by Majorana phases, which occurs even for two generations of neutrinos [7], CP violation in neutrino 
oscillations is a genuine three (or more) flavor effect, so it can be observed only when there is an 
interference between flavor oscillations involving at least two different Am^ and three mixing angles, as 
we will see below. We observe that the accelerator based neutrino oscillation experiments will provide 
the most promising opportunities to observe such CP violation. This will be extensively discussed in 
what follows. 

From the expression in Eq. (64), it is clear that there is no CP (or T) violation for three generation 
neutrino mixing in vacuum if 5 = or tt and 

• one (or more) of the mixing angles is zero, 

• two or more of the masses are degenerate. 

Moreover, note that it is impossible to observe CP violation in the disappearance channels {a = (3) 
since — ^ I'a is related to Ua '^a by CPT. Finally, there is no CP violation if the solar oscillations 
are averaged out. 

^Strictly speaking, even in vacuum, the oscillation probability for normal and inverted hierarchy differ and give a tiny 
effect, but we do not consider such a effect in this review, see [162]. 
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While AP defined in Eq.(64) vanishes for the disappearance channels, 

= 1 - 4\Ual\''\Ua2\^ Sin^ A21 - 4| [/c,2 n t/a3 1 ' sin^ A32 - 4 1 t/a3 H Ual \ ' sin^ A31 , (66) 

by CPT (hence no direct CP violation), the i/^ —>■ and Vj- — > v^. disappearance probabilities depends 
on cos 5, using our parametrization of the mixing matrix^. Given that if one knows cos 5 one also knows 
the magnitude of sin 5 up to a sign, can this be used to say anything about leptonic CP violation? 
In principle, the answer is yes. However it requires precision measurements in two out of the three 
disappearance channels so that one can determine precisely the magnitude of four independent elements 
of the mixing matrix. 

Suppose that one has measured \Ue2\1 jf^esl, |f^^j2| aiid|[/^3| precisely, then from unitarity one knows 
the moduli of all the elements of the mixing matrix. Construct the triangle which has sides with length 
|f^ei| I f^e2 1 1 f^^t2 1 and I ^/e3 1 1 ^/i3 1 ; ^s showu in Fig. 20. This triangle must close due to the following 
unitarity relation 

U:,U^, + Ul^U^2 + C/:3C/m3 = 0. (67) 

Now, twice the area of this triangle is the absolute value of the CP- invariant factor ^. Using the rule 
relating the length of the sides to the area of the triangle, we have 

J2 = ^[\Uel\\U^l\ + \Ue2\\U^2\ + \Ue^\\U^^\] X [-\U,,\\U^i\ + \U,2\\U^2\ + | f/e3 1 1 [^^3 1 ] 

A\Uel\\U^l\ - \U,2\\U^2\ + |^e3||^M3|] X [|^el||^Mll + \Ue2\\U^2\ ' |^e3||^M3|]- (68) 

The maximum possible difference in the lengths of the two sides, |C/ei||t^;ii| and |C/e2||t^//2| is |C^e3||t^/:i3|- 
At these extrema, the triangle fits on a line and has zero area and CP is conserved. Thus, to claim 
CP violation one must show at the required confidence level that ||C/ei||t^/ii| — |t42||C^/i2|| < |t^e3||t^/i3|- 
Given that |C/e3| is known to be small [12] and is yet unobserved, this will be a formidable challenge. 
The sign of J or equivalently the sign of CP violation cannot be determined from vacuum disappearance 
measurements. 

Note that the center of the circle in the triangle diagram bi-sects the | t/e3 1 1 ^/i3 1 side in the ratio of 
|t^e2p : \Uei\'^ from left to right. (This is equivalent to holding the "(1,2)" vertex fixed and rotating the 
"3" side about the same point.) Using the PDG parametrization of the lepton mixing matrix, the base 
of this triangle is S13C13S2 and the height C13S12C12C23I sm5\. The angle u) equals 5 or 27r — 5. 

The CP-invariant factor for any density matter is given by J{N) = (si2Ci2S23C23Si3C^3 sin 5)7v, where 
the mixing angles and CP phase arc their values in matter, obtained by writing U{N) in the form of 
Eq. (53). The matter value of the CP-invariant factor is related to the vacuum value as follows [163]: 

j( N) = Ami, A^ii J .g9^ 

'^''^ Ami, (N) Ami2 (N) Aml{N) ^^^^ 



•^Since |?7^ip, |?7^2p, iC^riP and |f/^2p depend on cos ^, e.g. |f/;:i2|^ = 0^3^12 + Si3«23«i2 - 2C12S12C23S23S13C13 cos5. 
^AU other unitarity triangles have the same area. 
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• — 




\Ue3\\Uf,s\ 

Figure 20: The unitarity triangle using the first and second rows of the mixing matrix. The lengths 
of each side are as labeled and twice the area of this triangle is the absolute value of the CP-invariant 
factor, I J|. The |f/ei| |f//ii| and | t/e2 1 1 f^^j2 1 vertex moves in a circle as the CP violating phase is changed. 

This identity guarantees that the AP in Eq.(64) is the same in matter as in vacuum for distances 
smaller than any of the matter or vacuum oscillation lengths. 

6.2 The oscillation probability — > z/g 

In this review, we mainly focus on the oscillation channel between electron and muon neutrinos because 
it is easier to create and detect these neutrinos compared to tau neutrinos. The drawing below shows 
schematically the relation among four possible channels. 

CP 

l^fi l^fi l^e 

T t $ T 

I'e ^fi I'e J^fi 

CP 

The horizontal (vertical) processes are related by CP (T) whereas the processes across the diagonals 
are related by CPT. The first row will be explored in very powerful conventional beams, Superbeams, 
whereas the second row could be explored in Neutrino Factories or Beta Beams. 
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6.2.1 Vacuum: — >■ i/g 

Let us first consider oscillation in vacuum ignoring matter effect. The transition probability v^, 
can be simply written as the square of a sum of three amplitudes, one associated with each neutrino 
mass eigenstate, as follows [164], 

= |2C/;3[/e3sinA3ie-*^^^ + 2[/;2C/e2sinA2i|^ (70) 

where the unitarity of the lepton mixing matrix [1] has been used to eliminate the U*iUei term and 
Ajk is used as a shorthand for the the kinematic phase, Am?-^L / . It is convenient to rewrite this 
expression as follows 

= i'atm + 2V^V^COs(A32 + (5)+Psol, (71) 

where as the notation suggests the amplitude \/ Patm depends only on Am|^ and yJP^oi depends only on 
Am2i. For propagation in the vacuum, these amplitudes are simply given by 

V-Patm = sin ^23 siu 2^i3 siu A31 

\fP^\ = COS ^23 COS 6^13 sin 26'i2 sin A21 ^ cos 6'23 cos ^13 sin 26112 A21 , (72) 

where ^13 and A21 are assumed to be small. In the amplitude V-fsob terms proportional to sin ^13 sin A2ie" 
have been neglected since they are of second order in the small quantities sin ^13 and A21. 

For anti-neutrinos 5 —5. Thus the phase between V-Patm and V-Psoi changes from (A32 + 5) to 
(A32 — 5). This changes the interference term from 

2V%^V^cos(A32 + (5) ^ 2v^V^cos(A32-5). (73) 

Expanding cos(A32 ± 5), one has a CP conserving part, 

2v^^^\/^cosA32Cos5, (74) 

and the CP violating part, 

=F2 a/ Patm \/P^\ sin A32 sin 5, (75) 

where - (+) sign is for neutrino (anti-neutrino). This implies that AP given in Eq. (64) can be rewritten 
as 

AP^p = P(i/^ ^ i/e) - P{V^ Pe) = -4:y/P^y/P~ism A32 sin S. (76) 

Note that this expression coincides with the exact formula in Eq. (64) because terms we neglected in 
Eqs. (71) and (72) do not contribute to the CP violating term. 
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Figure 21: The probability -P(z^^ Ue)'- (a) The components Patm and Psoi, (b) The full probability 
including the interference term for 5 = | and solid and dashed respectively. Since ^ = — | + 27r 
we have P{v^ — > z/g, 5 = ^) = P{u^ — > z/g, 5 = f ) so that these two probabilities can be considered to 
be -P(t'^ ^ Ue,5 = |) (solid) and P{p^ ^ = |) (dashed). The difference between these to curves 
demonstrates CP violation for this process. 

Therefore CP violation is maximum when A32 = {2n + l)f and grows with n since sin A21 grows 
with n. Notice also that, as expected, for CP violating term to be non-zero the kinematical phase 
A32 cannot be ravr. This is the neutrino counter part to the non-zero strong phase requirement for CP 
violation in the quark sector. 

Fig. 21 shows the components of ^(t'^ — * t'e) as well as the full probability for selected values of the 
CP phase 6 for neutrino energy E = 1 GeV. Unless otherwise stated, we fix the absolute value of Am"^^ 
to be 2.5x10-3 eV^, Amj-^ = 8.0 x 10"^ eV^, sin^^ia = 0.31 [12]. Since we can use CP and T to relate 
all the processes discussed earlier in section we have 



Pil^l, l^e, 5) = Piiy^ -> Z/e, -6) = P{Ue U^, -6) = P{Ue ^ t'/i, (5), 



(77) 



all given from Eq. (71). 

In the left panel of Fig. 22 we show APj^y in vacuum as a function of L for 5 = 0, 7r/2, tt and 3tt/2 
for sin^ 26'i3 = 0.05. As mentioned before, peak of |APi,p| occurs at A32 = (2n + 1)| growing linearly 
as n or L increases. 

For a fixed value of 6, as 613 becomes smaller, the absolute value of AP become smaller but this does 
not means that CP violation gets smaller. Consider the neutrino-anti-neutrino asymmetry, defined as 



A 



P - P 2v/Patmy^sin A32 sin S 

P + P Patm + 2 V-PatmV^ COS A32 COS 5 + Pgol ' 
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Figure 22: Examples of APj,p = P(z/^ 



(left panel) and the asymmetry AP/[P(fp 



z/e) + P(z^^ —>■ z^e)] (right panel) in vacuum as a function of distance for fixed value of energy, E = 1 
GeV and sin2 20i3 = 0.05. 



as a measure of CP violation. Then, at the vacuum oscillation maxima, A31 = {2n + l)f 



2 ' 



A 



atm 



P 



sol 



which is maximized when \^P^ = y/P^i- This occurs at 

sin^ 26*13 = cot^ 623 sin^ 26'i2 sin^ A21 ^ cot^ 6123 sin^ 20 



12 



Am|_ 



A^ 



(79) 



(80) 



At the first vacuum oscillation maximum, A31 = ^, the peak in the asymmetry occurs when sin^ 26'i3 ^ 
0.002. In the right panel of Fig. 22 we show the asymmetry in vacuum as a function of L for fixed 
value of sin^26'i3 = 0.05 and energy E = 1 GeV. As function of sin^ 26*13 the neutrino- anti-neutrino 
asymmetry is shown in Fig. 23. At the second vacuum oscillation maximum the asymmetry is maximum 
when sin^ 2^13 is 9 times larger. 

Another feature of the probability that is important is the existence of zero mimicking solutions [165]. 
When, 

-2/^cos(A32±5) (81) 

we have 



sol- 



(82) 



Thus, even if sin^2^i3 7^ 0, when this condition is satisfied non-zero 6*13 is impossible to observe. See 
Fig. 23 for examples of these "zero mimicking solutions" . The maximum value of sin^ 2^13 for which 
there are zero mimicking solutions is 4 times the value given in Eq. (80). Luckily, the zero mimicking 
solutions appear at different values of 6 if E/L is changed or by switching to anti- neutrinos. 
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Figure 23: The left panel (a) shows the asymmetry between neutrino and anti-neutrinos as a function 
of sin^26'i3. The peak occurs when VPatm = V-Psoi at sin^2^^i3 ^ 0.002. The right panel (b) shows the 
location of the zero mimicking solutions at the first oscillation maximum. 

6.2.2 Matter: 

In describing neutrino propagation in the Earth it is very important to take into account the effect of 
matter [18, 143, 19] because it can induce a fake CP violating effect even if the CP phase is zero or vr. 
In matter of constant density, the two component appearance probability, P^p, using the mixing angle 
and mass squared difference in matter, 6n and Am^, is given by 

P^pp = sm^2er,sm^A^ (83) 

where A^v = Amj^L/AE. Since Am1fSm29N is independent of the density of matter i.e. it is an 
invariant, then 

i'.^ = sin^2^o(^^)A^. (84) 

Note, that this expression does not depend on the mixing angle in matter, 6^ and that Am^ is given 
by 

Am% = ^J {Ami cos 200 - 2y/2GFN^Ef + {Ami sin 2^o)^ (85) 

where Gi;' is the Fermi constant and N^. is the number density of electrons. Except near the resonance, 
Am% = Ami - '^V^GpNeE is a good approximation. For matter effects to significantly alter the 
appearance probability the following condition must be satisfied, with the left and right hand side 
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significantly different 



sin^ A 



N 



sin^ Ac 

A2 
^0 



(86) 



That is, matter must significantly alter the Am^ so that Amj^ ^ AttIq and the baseline of the experiment 
must be a significant fraction of the oscillation length in matter or vacuum whichever is shorter so that 
either sin Ajy 7^ A^v and/or sin Aq 7^ Aq. 

For the three neutrino case, nature has chosen two small parameters, sin^ ^13 < 0.04 and Am|]^/ Aml^ r 
0.03 [12]; this allows us to factorize the three neutrino case into a product of two neutrino cases and 
therefore the individual Am^ in matter become 



Am^ 



Am^ 
Am^ 



■31 1 

2 I 
'21 1 

2 I 
■32 I 



2V2GfN^E 



(87) 



Am32. 



In Fig. 25 we have plotted the exact mass squared differences in matter and the approximation given 
in Eq. (87) and indeed the approximation is a good one. 
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Figure 24: Examples of AP = P(z^e ^p) ~ P{i^e ^p) (left panel) and the asymmetry AP/[P{Uf, — >■ 
z/^j) +P(z/e z/^J] (right panel) in matter as a function of distance for fixed value of energy, E = 1 GeV 
and sin^2^^i3 = 0.05. For simplicity, we assume a constant electron number density Ne = 1.5 mol/cc. 



Thus the ^/Patm and ^/Psoi in matter are simply given by 

\/ -Patm 



sin(A3i - aL) 
sm6'23 sm26'i3 — ^r— A31, 



sol 



(A31 - aL) 

„ . ^„ sin(aL) ^ 
cos 6^23 sm26'i2 , A21, 
[aL) 



(88) 
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Figure 25: The neutrino mass squared differences in matter as a function of energy (negative energy is 
anti- neutrinos). The density of matter is 3 g-cm~^ and vn\ at zero energy (vacuum) is arbitrarily chosen 
to be eV^. 



where a = GFNe/V2 which is approximately (3500 km)~^ for pYe = 3.0 g ■ cm~^. The relative phase 
(A32 + 6) between ^/Patm and y/Psoi in Eq- (71) remains unchanged. It is clear from \/P^ that relative 
size between the kinematic phase A31 compared to aL determines the effects of matter provided at least 
one of these is bigger than 7r/4. Since L is common to both A31 and aL provided it is a significant 
fraction of an oscillation (either in matter or vacuum) then the comparison is between Amli/E and 
a = G pN This implies that the larger the energy of the neutrino, the larger the matter effect 
provided the baseline is the same fraction of an oscillation length. 

In Fig. 26 we show the effects of matter on Patm and Psoi for fixed energy and varying the baseline 
and for fixed baseline and varying the energy. For fixed energy, one sees that the amplitude of the 
matter effect does not increase with distance but the shift in the peaks increases in proportion to the 
baseline. For fixed baseline, the amplitude of the matter effect gets smaller with lower energy, also 
proportionally. Thus the biggest matter effect is at the first oscillation peak. 

Thus in matter to leading order 



sin^ ^23 sin^ 2^13 



sin^(A3i - aV) 



(A31- 
+ sin 2^23 sin 2^13 sin 2^^i2 



sin(A3i 



^31 



aV) 



(A31 - aL) 



sin(aL) ^ , . ^, 

^31 , \. A21 cos A31 + S) 

[aL) 



+ COS^6'23 



9 ^ sin^(aL) , , 



(89) 



where the first (last) term is the atmospheric (solar) probabilities and the middle term is the interference 
between the atmospheric and solar contributions. The phase of the interference has been written as 
( A31 + (5) instead of ( A32 + 5) since the difference between these two is small. Also, sin 2^13 is often 
written as 2 sin ^13 or even 2^^i3, again the difference is higher order. We have allowed for A31 and A32 
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Figure 26: The components of the probabihty P{i'i_i Ve) in matter: The left column gives Patm and 
Psoi, whereas the right column gives the full oscillation probability including also the interference term 
between the atmospheric and solar amplitudes. The first row is for fixed energy, 0.6 GeV, and varying 
the baseline, whereas the second row is for fixed baseline, 1200 km, and varying the energy. 

to be positive or negative in this derivation, thus this expression is valid for both normal, Amg^^ > 
and inverted, Amg^ < 0, hierarchies. 

For the CP conjugate process, z/^ Ve-, in matter both a —a and S ^ —S is required. Also 
for the T conjugate process, u^. —>■ z/^, replacing L with —L including in the kinematic phases. A, is 
necessary. Note also that the CPT conjugate process, Ue —>■ i^fi, is obtained by changing the sign of a, 5 
and L. If one changes the sign of a again, since this is equivalent to having propagation in anti-matter, 
the probability is once again the same as that for z/^ as it must from CPT invariance, since 

P{^fj. ^ i^e) I matter = Pi^^e ~^ '^/i ) | anti-matter • (90) 

This CPT relationship is true for any symmetric matter profile, a{x) = a{L — x), where x is the position 
along the neutrino path and L is the baseline of the experiment. 
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In fact there are a number of relationships between these probabihties [166, 167]. Here we discuss 
some, especially those that hold for the case of inverting the hierarchy. The relationship between CP 
conjugate pairs is given by 

^ P{iye ^ i^fi;-Aml^,+Aml^,TT + 5,a). (91) 

The CPT relationship including the effects of switching the hierarchy is 

P {i^fi ^ ^e] ^1^31, Amli, 5, a) = P^Ue ^ i^^i] —Aml^, ~Aml^,2TT — 5, a) 

fa P(z/e ^ z/^; -Amg^, +Am2i, TT - 5, a). (92) 

This relationship is demonstrated in Fig. 27 by the dotted horizontal lines. In fact, there are many 
such relationships, e.g., 

P{u^ ^ Ue] Amg]^, Am2i,5, a) = P{u^ ^ Ue', —Am^^, —Am2i,2n — 5, —a) 
= P{u^ I'e] —Amli, —Amli, 5, a) =P{iJ^^Ue] Am\i, Am\i,2Tx — 5, —a). (93) 

The easiest way to derive such relationships is to look at the evolution equation, Eqs. (56) and (59), 
and its complex conjugate for neutrinos and anti- neutrinos using a symmetric matter profile and noting 
that U*{2ti-5) = U{6). 



V, ^ v„ E,=13 GeV L=3000 km 




<P(i^)> % 



Figure 27: P{u^ ^ Ve) versus P{u^ v^) and Piv^ — ^ Ve) versus P{ve ^^i) in matter. In vacuum 
the four ellipse overlap, one for CP conjugate pair plus one for T conjugate pair times two, one for the 
normal (+) and one for the inverted (— ) hierarchies. When matter effects are turned on the four ellipses 
split to form a baseball diamond with the vacuum ellipse is the pitcher's mound. All the parameters 
are held fixed except the CP- violating phase, 6 and sin^ 26*13 = 0.05. The dotted lines demonstrate the 
approximate relationship of Eq. (92). 
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We show in Fig. 24 the same quantities shown in Fig. 26 but in the presence of the matter effect. 
As we can see, the matter effect can make AP^^p as well as the asymmetry significantly different from 
zero even if the CP phase is zero or tt. 

When the matter effect is significant the quantity defined in Eq. (76) is no longer a good measure of 
the intrinsic CP violation, since the matter effect makes this quantity significantly different from zero 
even if 5 = or vr. Therefore we define 

5Pcp = AP,p(5 = 7r/2)-AP,p(5 = 0), (94) 

which reduces to Eq. (76) in vacuum. We note that in this quantity the matter effect is approximately 
canceled so that only the intrinsic CP violation effect remains, to a good approximation. 




E/L [GeV/km] 

Figure 28: Iso-contours of 5Pcp = AP^p{6 = 7r/2) - AP^p(5 = 0) for sin^ ^13 = 0.05 in the E/L - L 
plane in matter for normal hierarchy. 

In Fig. 28 for a given baseline distance the matter density was assumed constant and estimated 
from the Preliminary Reference Earth Model (PREM) [168] by taking the average along the neutrino 
trajectory. The baseline and energy corresponding to T2K and NOz/A are indicated by the solid circle 
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and square, respectively. The positions of the peaks of the 5 "peninsula" regions correspond to, from 
right to left, 1st, 2nd, 3rd, 4th and 5th oscillation maxima. The best places to search for CP violation 
effects in the probabilities are those within the regions with higher 5Pcp values. Since the neutrino flux 
decreases as as the baseline gets larger one should, in practice, try to measure the first or first few 
oscillation maxima. We note that in the region with L ~ 7000 — 8000 km no significant CP violation is 
expected: at such "magic" basehnes [169] the dependence of the probability on the CP S and the solar 
oscillation parameters disappears (see [170] for its physical interpretation). 

7 Near Term Experiments 

Here we discuss the prospects of observing CP violation by the following two near term experiments^, 
T2K and NOuA: 

• T2K [171]: Initially this experiment will send a beam of z/^'s from the JPARC facility to the 
existing 50 kton Supcr-Kamiokande (SK) detector. The beam will be aimed 2.5 degrees off the 
JPARC to SK axis resulting in a mean neutrino energy of 0.63 GeV with an approximate Gaussian 
spread of 0.13 GeV at SK. The source detector distance is 295 km with a mean matter density 
of 2.3 g cm^'^. The beam power will be ramped over time to 0.75 MW. Only neutrino running is 
expect during the initial phase, Phase I. In the future, there is a possibility to upgrade the beam 
power to as high as 4 MW and to include anti-neutrino running. 

• NOz/A[172, 173]: The NOz/A detector will be placed 12km off-axis of the existing NUMI neutrino 
beam line at a distance of 810 km. The mean neutrino energy in this configuration is 2.0 GeV with 
an approximate Gaussian spread of 0.30 GeV. By the time the 25kton liquid scintillator detector 
has been completed the NUMI beam power is expected to be greater than or equal to 0.4 MW. 
This allows for the running in both neutrino or anti-neutrino mode. The mean matter density for 
this configuration is 2.8 g cm~^. In the future, there is a possibility to upgrade the beam power 
up to as high as 2 MW. 

Given the current best fit value for Am^tm — 2.5 x 10~^ eV^ from the global analysis [174] both of 
these experiments are at a mean energy which is above the first vacuum oscillation maximum (VOM) 
energy. Also, because of the longer path length for the NOz/A experiment the matter effects are three 
times larger for NOz/A than the T2K experiment. 

*At the time of writing this review, the beam hne for T2K is under construction with first beam to the existing 

Super-Kamiokande detector is expected in 2009. The NuMI beam line is operating for the MINOS experiment but the 
NOz^A detector has still to get final approval from DOE to begin construction and is expected to come on line after T2K 
by at least one year. 
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7.1 Dependence of probabilities on ^13 and 6 



First let us look at the how the oscillation probabilities as well as CP violation (the difference of 
probabilities) for these two experiments depends on 613 and 6. Fig. 29 shows the iso-probability contours 
for the T2K experiment for both hierarchies in neutrino and anti-neutrino running. The expected 
sensitivity for the first phase of running of T2K is approximately equal to the 0.5% contour in neutrino 
running. There is a small difference between the contours for the two hierarchies due to matter effects 
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Figure 29: Iso-contours of P{jyfj, — ^ J^e) and -P(t'^ z^e) in matter for T2K experiment. The solid 
(dashed) curves correspond to the normal (inverted) hierarchy. In Phase I, T2K is expected to reach a 
sensitivity which is close to the 0.5% contour for neutrinos. 



and the fact that the mean neutrino energy is above vacuum oscillation maximum. Since the beam 
and fake backgrounds contribute approximately 1%, further improvement on the sensitivity will be 
challenging. No anti-neutrino running is expected in the first phase since the event rate is appreciable 
smaller than in neutrino running. Of course, anti-neutrino will be included as the beam power is ramped 
up above the initial phase. Due to the increased background and smaller statistics, compared to neutrino 
running, a reasonable estimate of the sensitivity in anti-neutrino running is the 1% contour of the anti- 
neutrino plot in Fig. 29. Fig. 30 shows the iso-probability contours for the NOi/A experiment for both 
hierarchies in neutrino and anti-neutrino running. The expected sensitivity for neutrino running is also 
approximately 0.5% and approximately 1% for anti-neutrino running. Again, further improvements 
will be challenging due to backgrounds. Notice the sizable difference in the sensitivity between the 
two hierarchies [175], primarily due to the matter effects. The combined running of neutrinos and 
anti-neutrinos will make the combined sensitivity much less dependent on the CP violating phase, 6. 

In order to have some idea about the magnitude of the intrinsic CP violation (coming from 6) 
which can potentially be observed in these two experiments, we show in Fig. 31, the iso-contours of 
APj^p = Pi^i'^ I'e) — Pipfi ~^ ^e) in the plane of 5 and sin^ 2^13 by switching off the matter effect. 
Note, that this quantity is just the CP-invariant factor multiplied by a constant which depends on the 
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Figure 30: Iso-contours of P{iy^ — > I'e) and Pii^fi Ve) in matter for NOz/A experiment. The solid 
(dashed) curves correspond to the normal (inverted) hierarchy. The NOi^A experiment in its initial 
phase is expected to reach a sensitivity close to the 0.5% contour in neutrinos and approximately the 
1% contour in anti-neutrinos. 



neutrino energy and baseline, see Eq. (64). 

Note also that in vacuum this quantity does not depend on mass hierarchy so that the contours 
are symmetric with respect to the line 5 = vr except for the sign difference. The positions of the CP 
conserving solutions in the (sin^ 2^13, sin 5) plane are shown in Fig. 32. 

However, this quantity is not directly observable as the actual positions of the APj^p contours are 
modified by matter effects, as shown in Fig. 33 for T2K and Fig. 34 for NOz/A. As expected, the 
modifications due to the matter effect, which can be seen by the asymmetry of contours with respect 
the line 5 = vr, are smaller for T2K than for NOz/A. 

In Fig. 32 we have show the location of the CP conserving solutions of one hierarchy in the 
(sin^ 26'i3, sin 5) plane of the other hierarchy. These locations depend on the matter effect for the 
experiment and thus differ for T2K and NOz/A. In order to claim the observation of CP violation, one 
must be sufficiently away (with a required CL) not only from the line sin 5 = but also from the 
dot-dashed (T2K) and dashed (NOz/A) curves shown in these plots unless the mass hierarchy is known. 

For baseline L < 1500 km or so, one can expand APi,p(a) in matter as a Taylor series about the 
vacuum AP,^p(0) as follows, 

AP,p(a) = AP,p(0) ±4(aL)G(A32)[Patm(0) + VPatm(0)Psol cos A32 COS 5], (95) 

where G{x) = 1/x — cot a; and the + (-) sign is for the normal (inverted) hierarchies. This simple 
expression gives the relationship of Figs. 33 and 34 to Fig. 31. 
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Iso-Contors of AP - [%] in Vacuum (L=295 km, E = 0.65 GeV) Iso-Contors of AP - [%] in Vacuum (L=810 km, E = 2.0 GeV) 




5/7C 6/7C 

Figure 31: Iso-contours of AP^p = P^v^ u^) — P{j^^ v^.) in the 5 and sin^ 2^13 in vacuum for T2K 

(left panel) and NOi/A (right panel). 

7.2 Bi-probability plots 

Of course, the neutrino and anti-neutrino probabilities are correlated. This correlation depends on the 
CP violating phase, 5, and the mass hierarchy. Therefore it is useful to use bi-probability diagrams. In 
the bi-probability space spanned by P{vi_i — > Ve) and -P(t'^ i^e), for a given ^13 (fixing all the mixing 
parameters except for 6), energy and baseline, the variation of 6 from to 2tt gives a closed trajectory 
which is an ellipse^. In Figs. 35 and 36, we give the bi-probability diagrams for both T2K and NOz/A, 
respectively. For T2K the separation between the hierarchies is small but non-negligible whereas for 
NOi/A the separation is much larger with complete separation occurring for sin^ 2^13 = 0.11. The value 
of ^13 for which the two hierarchies separate is called the critical value, For T2K and NOi/A the 
bi-probability diagrams are shown for two energies, one such that A31 ^ tt /2 which is known as vacuum 
oscillation maximum (VOM) and at an energy above VOM. At VOM the ellipses in the bi-probability 
plot are squashed to lines as the probabilities become independent of cos 5 since COSA31 = 0. 

7.3 Parameter Degeneracy 

Unfortunately the measurement of Pii^fj. —>■ Ve) and P{i'^ —>■ Ve) does not determine the values of ^13 
and 5 uniquely due to the possibility of drawing more than one ellipse through any (P, P) point. In 
general four such ellipses can be drawn assuming sin^ 622, is known uniquely. Eight if only sin^ 2^23(7^ 1) 
is known. Fig. 37 demonstrate this point. 

The cross in the left hand panel of these two figure, are the transition probabilities for the normal 
hierarchy with sin^26'i3 = 0.05 and 5 = 7r/4 for both T2K and NOz/A. The right hand panels show 
the four allowed solutions in the sin^ 26'i3 v sin 5 plane that are consistent with this point in the (P, P) 

^We note that the trajectory is exactly elliptic as it was shown [176] that even in matter, oscillation probabilities take 
the form oi P — A cos 5 + B sin 5 + C where A, B and C are some constants which depend on mixing parameters as well 
as experimental parameters such as energy and baseline. 
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Figure 32: The positions of the CP conserving solutions in the (sin^ 2^13, sin 6) plane: (a) for the Normal 
Hierarchy, the curves above sin 5 = for T2K (dot-dashed) and NOz/A (dashed) are the locations of the 
CP conserving solutions for the Inverted Hierarchy, (b) is the equivalent plot for the Inverted Hierarchy. 
Without determining the hierarchy, a given experiment cannot claim the observation of CP violation 
unless it can exclude both of lines with its label at the required confidence level in one of these plots. 



plane. Clearly, the four solutions are related to one another: 



• Within the same hierarchy, the intrinsic degeneracy [177]. 
The two solutions have approximately the same sin 5, but different signs for cos 5 i.e. 5 and n — 5. But 
the values of 613 differ by 

A6'i3 = cos 6 sin 26'i2 A21 cot A31 



Iso-Contors of AP - [%] in Matter for NH (L=295 km, E = 0.65 GeV) Iso-Contors of AP - [%] in Matter for IH (L=295 km, E = 0.65 GeV) 




Figure 33: Iso-contours of AP^p = P{v^ — > z/g) — P^v^ —* Ve) in the 5 and sin^ 2^13 for T2K for the 
normal (left panel) and the inverted (right panel) mass hierarchy, with the matter effect included. 
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Iso-Contors of AP - [%] in Matter for NH (L=810 km, E = 2.0 GeV) Iso-Contors of AP - [%] in Matter for IH (L=810 km, E = 2.0 GeV) 
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Figure 34: Iso-contours of AP^p = Pi^u^ v^) — P{.^n ^ ^e) in the 5 and sin^ 26'i3 for NOz/A for the 
normal (left panel) and the inverted (right panel) mass hierarchy, with the matter effect included. 
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Figure 35: The bi-probability P{i'f_i —>■ Ve) diagram for T2K including matter effects. 
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Figure 36: The bi-probability -P(i^^ Ve) diagram for NOz/A including matter effects. 
At vacuum oscillation maximum, A31 = 7r/2, these two different values of 6^13 coincide, see Fig. 37. 



• Between the two hierarchies[175]. 
The allowed values for sin^ 2^13 are the same for both hierarchies but the sin 5's differ between the 
normal (NH) and inverted (IH) hierarchies as follows, sin5|7VH — sin5|/// = 26'i3/6'j3*, see [178]. Where 
is the largest value of for which the two regions in the bi-probability diagram overlap^^, see 
Fig. 37. For T2K and NOz/A this corresponds to 



sin(5|7VH — sin5|/// = < 



0.47 



1.42 



sin^ 26'i3 
0.05 



sin^ 2^13 
0.05 



T2K 



NOz/A 



(97) 



The coefficient is 3 times bigger for NOz/A than T2K primarily because the baseline for NOz/A is 
approximately 3 times that of T2K so that the matter effects are 3 times larger. Notice that the matter 
effects are smaller but not negligible for T2K, see Fig. 32. Clearly the combination of T2K and NOz/A 
can be used to determine the hierarchy at least for large values of ^13. 

If sin^26'23 7^ 1 then there is a further degeneracy [179] associated with whether 



sin^ 6 



1 T VI - sin" 2^23 



23 



The value of f^^s for any experiment is given by 



9crit 
13 



TT^ sin26'i2 Amji 
8 tan6'23 ^7713^ \y 1 — A31 cot A31 



/{aL). 



(98) 



(96) 



55 
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Figure 37: (a) The bi-probability plot for T2K (top) NOz^A (bottom), showing the allowed regions for 
the normal (dashed) and inverted (dot-dashed) hierarchies. The cross is the assumed measured value 
of {P,P). The value of Oi^ for which the two hierarchies separate is labeled as the critical value, ^J"*. 
(b) The four solutions for the assumed point in the sin^26'i3 vs. sin 5 plane. Since T2K is near vacuum 
oscillation maximum all four solutions have approximately the same value of sin^ 2^13. But the solutions 
for the normal and inverted hierarchies occur at different values of sin 5. Since NOz/A is above vacuum 
oscillation maximum, there are two distinct allowed values of sin^ 2^^i3. Again the solutions for the same 
hierarchy have the same value of sin 6. But the solutions for the normal and inverted hierarchies occur 
at different values of sin 6. 
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The simultaneous presence of these three independent 2-fold degeneracies leads up to 8-fold degen- 
eracies [180]. Resolving all three of these degeneracies will be challenging but, except for the hierarchy 
degeneracy, they do not affect much whether or not one can claim the observation of leptonic CP 
violation. 

For the mass hierarchy degeneracy this ambiguity can affect whether or not one can claim CP 
violation especially for larger values of ^13. At larger values of ^13, NOi/A alone may be able to 
determine the hierarchy if sign(Am|i) sin 5 is close to one. Otherwise a combination of T2K and NOi/A 
is our best option for determining the hierarchy with only these two experiments. Even if the hierarchy 
is not determined there is still a substantial region of (^13, 5) space where not knowing the hierarchy does 
not preclude being able to claim the observation of leptonic CP violation. However there is a region in 
both experiments where one hierarchy is CP violating and the other CP conserving, fortunately these 
troublesome regions occur in different locations for the two experiments. This makes the combination 
of T2K and NOi/A a very powerful tool to observe leptonic for large ^13. 

7.4 Beyond the First Oscillation Maximum: 

To untangle all of the degeneracy associated with the neutrino mass hierarchy, the quadrant of $23 
and the sign of cos^cp discussed in the previous subsection, it is probable that experiments beyond 
the first atmospheric oscillation maximum will be necessary. First, let us discuss what happens in 
vacuum [181] and then add matter effects. In vacuum, the atmospheric amplitude given by Eq. (72) 
has the same magnitude at each of the successive maxima whereas the solar amplitude at the n-th peak 
is {2n — 1) times larger than at the first peak as long as A21 is less than 1. Thus, the magnitude of 
the ratio of the solar to atmospheric amplitude grows as (2n — 1). Hence, the relative magnitude of the 
probabilities associated with the atmospheric and solar Sm^ and their interference in the total i/^ — > i/e 
transition probability, the three terms in Eq. (71), changes at the successive peaks. At the n-th peak, 
the interference term, and hence CP violation, grows as (2n — 1) compared to the first peak whereas the 
term associated with purely the solar Sm^ grows as (2n — 1)^. This change in the relative magnitude of 
the three terms can be exploited to untangle the degeneracies. 

In vacuum, the transition probability z/^ — > is the same if one goes to the successive peaks by 
lowering the energy or increasing the baseline or some combination of the two (as long as E/L is kept 
to be the same). However, as we will see shortly, in matter, the — > transition probability depends 
on exactly how one goes to the successive peaks. E.g. the transition probability at the second peak is 
different for lowering the energy by a factor of three at the same baseline than increasing the baseline 
by a factor of three at the same energy. To understand this wc need to understand how matter effects 
change the atmospheric and solar amplitudes which is most easily done for a constant matter density. 
The solar amplitude changes little except at very large distances whereas the atmospheric amplitude 
is substantially modified by matter effects, see Eq. (88). The change in the atmospheric amplitude 
depends on whether one goes to the successive peaks by varying the baseline holding the energy fixed 
or varying the energy and holding the baseline fixed. 
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Figure 38: Examples of the bi-probability -P(t'^ — > t^e) diagram corresponding to near the second 
oscillation maximum including matter effects. 

For fixed energy and varying baseline, the modification of the atmospheric amplitude can be under- 
stood as just a change in the mixing angle and which depend on the energy and matter density 
but are independent of the baseline, see section 6.2.2. So the amplitude of the oscillation is the same at 
all successive peaks and the position of n-th peak is just {2n — 1) times the distance to the first peak. 
See top left panel of Fig. 26. 

For fixed baseline and varying energy, again one can think of the modification as just changing the 
mixing angle and Sm"^ but these parameters are energy dependent and approach the vacuum values 
as the energy gets smaller and smaller, see Fig. 25. Again, see the discussion in section 6.2.2. So the 
matter effects are reduced by (2n — 1) in amplitude at the n-th peak since the energy at the n-th peak 
is reduced by (2n — 1) compared to the first peak. Therefore, matter effects become increasingly less 
important at the successive peaks. This can be see in the bottom left panel of Fig. 26. 



7.5 Possible Experiments 

A possible extension of the T2K experiment, called T2KK [182, 183] (see also [184, 185]), consists of 
building two new large detectors, one at the first oscillation peak at Kamioka, Japan and the other in 
Korea, near the second oscillation peak. Both detectors would be at the same off-axis angle so they 
would see the same neutrino energy spectrum but the Korean detector would be three times further 
from the source than the Kamioka detector. This is an example of getting to the second peak by varying 
the baseline holding the neutrino energy fixed. In the left panel of Fig. 38 we show the examples of 
bi-probability plots corresponding to such an experimental set up, which should be compared with the 
right panel of Fig. 35. From these two plots, we confirm that the oscillation probabilities as well as its 
dependence on S are significantly enhanced at Korea due to the increase of the solar term. In Fig. 39 we 
show the iso-contour of APj^p = Piy^ Ve) — PipiJ. J^e) in the 5 — sin^ 2^13 plane for the normal (left 
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Iso-Contors of AP - [%] in Matter for NH (L=1050 km, E = 0.65 GeV) Iso-Contors of AP - [%] in Matter for IH (L=1050 km, E = 0.65 GeV) 




Figure 39: Iso-contours of AP^^p = ~^ ^e) — Pipn ^e) in the 5 and sin^2^^i3 for T2KK (L=1050 
km) for the normal (left panel) and the inverted (right panel) mass hierarchy, with the matter effect 
included. 



Iso-Contors of AP - [%] in Matter for NH (L=810 km, E = 0.6 GeV) Iso-Contors of AP - [%] in Matter for IH (L=810 km, E = 0.6 GeV) 




Figure 40: Iso-contours of APj,p = Pi^u^ Ue) — P{i^^ z^e) in the 6 and sin^26'i3 for the NOz/A 
baseline (L=810 km) but for the neutrino energy E = 0.6 GeV for the normal (left panel) and the 
inverted (right panel) mass hierarchy, with the matter effect included. 

panel) and the inverted (right panel) hierarchy, which should be compared with Fig. 33. We observe 
that the qualitative behaviour of these plots are similar but magnitude of AP^^p for larger baseline is 
significantly greater. By studying the neutrino spectrum at both detectors one could untangle all of 
the degeneracy [183]. 

Another possibility is to build a second large off-axis detector near the NOz/A detector such that the 
energy of the neutrino beam is one third of that seen by the NOz/A detector [186]. Both detectors would 
be approximately at the same distance so this is an example of getting to the second peak by varying 
the neutrino energy and holding the baseline fixed. At the second peak, CP violation is three times 
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larger than at the first peak and matter effects are three times smaller. Therefore, counting experiments 
at both peaks could resolve the entanglement of the mass hierarchy and CP- violation. 

In the right panel of Fig. 38 we show an example of bi-probability plot corresponding to this case 
which should be compared with Fig. 36. We confirm that the effect of the CP phase (i.e., the size of 
the ellipses) become larger and the matter effects get smaller. The corresponding iso-contour plots of 
AP,^p, shown in Fig. 40, also support this feature. 

A third possibility is to build a large detector on-axis at a distance of 1000 - 1500 km from the 
source and to measure the oscillated neutrino spectrum over a wide neutrino energy covering the first, 
second and possibly the third atmospheric oscillation peaks [187, 186]. This again is an example of a 
fixed baseline and varying the neutrino energy. The hierarchy would be determined by studying the 
first peak and CP- violation from studies of the second peak. 



As wc have seen, neutrino oscillations are insensitive to the absolute scale of neutrino masses. While 
the origin of the latter remains a theoretical mystery, one expects on general grounds that neutrinos 
are Majorana fermions, and that this accounts for their relative lightness with respect to the other 
fundamental fermions. 

The search for lepton number violating processes such as f3f3ou opens the way to probe the basic 
nature - Dirac or Majorana - of neutrinos, and also to probe the CP violation induced by so-called 
Majorana phases [7]. As already mentioned, the latter does not show up in conventional neutrino 
oscillation experiments [24, 25, 26] but can affect neutrinoless double beta decays and electromagnetic 
properties of neutrinos [23, 28, 32] [188, 189]. 

The most direct mechanism engendering PPou is the so-called "mass mechanism" involving the 
exchange of massive Majorana neutrinos. The associated amplitude is proportional to 



where the K^i are the first row in the lepton mixing matrix. An important feature of this amplitude 
is that it involves the lepton-number-violating propagator in Eq. 21 and, as a result, none of the Kei 
factors appear with complex conjugation, hence the effect of Majorana phases and the possibility of 
destructive interference among amplitudes arising from different neutrino types. 

Such destructive interferences may also take place without CP violation. This is what happens for 
example in the case of a pure Dirac neutrino, the exact cancellation coming in this case from the phase 
present in Eq. (15). Related concepts are those of partial cancellations due to symmetry. One example 
is the case of a Quasi-Dirac neutrino, made up of one active and one sterile neutrino, almost degenerate 
in mass due to the nearly exact conservation of standard lepton number [29]. Mass splittings between 
the two neutrinos may be radiatively calculable in some gauge models [30] and lead to active-sterile 
neutrino oscillations. Depending on the size of the mass splitting, these active-to-sterile conversions 
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may affect primordial Big Bang Nucleosynthesis [190, 191]. A different concept is that of a pseudo-Dirac 
neutrino made up of two active neutrinos, nearly degenerate due to the approximate conservation of 
some non-standard combination of lepton numbers [31]. 

The significance of neutrinoless double beta decay stems from the fact that, in a gauge theory, 
irrespective of which is the mechanism that induces P(3oi^ , it necessarily implies a Majorana neutrino 
mass [23], as illustrated in Fig. 41 This is specially relevant given the fact that gauge theories 
bring in other possible ways of inducing /3/3oi/ and it is conceivable that one of such mechanisms might, 
perhaps, give the leading contribution. Hence the importance of searching for neutrinoless double beta 
decay. Quantitative implications of the "black-box" argument are strongly model-dependent, but the 




V V 



Figure 41: Neutrinoless double beta decay and Majorana mass are theoretically equivalent [23]. 

theorem itself holds in any "natural" gauge theory. 

Now that oscillations are experimentally confirmed we know that there is a contribution to (3(3oi, 
involving the exchange of light Majorana neutrinos, the " mass- mechanism" . The corresponding ampli- 
tude is sensitive both to the absolute scale of neutrino mass as well as the two Majorana CP phases 
that characterize the minimal 3- neutrino mixing matrix [7], none of which can be probed in oscillations. 

Fig. 42 shows the estimated average mass parameter characterizing the neutrino exchange con- 
tribution to PPoiy versus the lightest neutrino mass. The calculation takes into account the current 
neutrino oscillation parameters in [12] and latest nuclear matrix elements of [192] and compares with 
experimental sensitivities. 

The most stringent lower bounds on the half-hfe of PPo„ decay were obtained in the Heidelberg- 
Moscow ^^Ge [193] and CUORICINO ^^^Te [194] experiments: 

Ti/2C^Ge) > 1.9 • 10^5 years, Tl'^/^C^^Te) > 1.8 • 10^^ years. (99) 

Using recently calculated nuclear matrix elements with significantly reduced theoretical uncertainties 
[192] from these data the following upper bounds for the effective Majorana mass m/j/j = {rrii,) can be 

^^Such black-box theorem does not exist for the case of lepton flavour violation, which may proceed in the absence of 
neutrino mass [80]. 
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[eV] m3 [eV] 

Figure 42: Neutrinoless double beta decay amplitude versus current oscillation data [192]. 



inferred 

\^/3f3\ ^ 0.34 eV (Heidelberg — Moscow) 

\mi3f3\ < 0.55 eV (CUORICINO). (100) 

The Heidelberg group, which includes a few authors of the Heidelberg-Moscow collaboration, recently 
claimed [199] evidence for the (3(3o^ decay of ^^Ge with T^J^ = (0-69 - 4.18) ■ 10^^ years at the 4.2a 
confidence level. Using the nuclear matrix element obtained in Ref. [192], from this data one finds for 
the effective Majorana mass the range 0.23 eV < \Tni3j3\ < 0.56 eV. Future experiments will extend the 
sensitivity of current /S/Sou searches and provide an independent check of this claim and even go further 
in sensitivity, as seen in the figure [201]. 

The left (right) panel in Fig. 42 corresponds to the cases of normal (inverted) neutrino mass spectra. 
In these plots the "diagonals" correspond to the case of quasi-degenerate neutrinos [70], which give the 
largest f3(3ou amplitude. In the normal hierarchy case there is in general no lower bound on the (3(3ou 
rate since there can be a destructive interference amongst the neutrino amplitudes. In contrast, the 
inverted neutrino mass hierarchy implies a "lower" bound for the jS/Soi, amplitude. 

A normal hierarchy model with a lower bound on /5/3oi/ is given in Ref. [132]. An interesting feature 
is that the lower bound obtained depends on the value of the Majorana violating phase 0i, as indicated 
in Fig. 43. Note that the lines in dark (red) and grey (green) of the left panel correspond to normal 
and inverse hierarchy, respectively. An alternative model based on the A4 flavour symmetry has been 
suggested [133] which implies a lower bound on the neutrinoless double beta decay rate, corresponding 
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Table 2: Sensitivities of future Oz//5/3-decay experiments to the effective Majorana neutrino mass cal- 
culated with the RQRPA nuclear matrix elements M^^{A,Z) of Ref. [192]. For the axial coupling 
constant qa the value qa = 1-25 was assumed. T^Jg"^^^ is the maximal half- life, which can be reached 
in the experiment and mpp = {niv) is the corresponding upper limit of the effective Majorana neutrino 
mass. 



Nucleus 


Experiment 


Source 


. 11 7 y - - 'r* T» r n 

^1/2 ' N 


Ref. 




\mpp\ [e 




GERDA(I) 


15 kg of ^"^''Ge 


3 10^5 


[200] 


3.92 


0.27 




GERDA(II) 


100 kg of ^""''Ge 


2 10^6 


[200] 


3.92 


0.10 




Majorana 


0.5 t of ^"''Ge 


4 10^7 


[195] 


3.92 


0.023 




SuperNEMO 


100 kg of ^"^^e 


2 10^6 


[196] 


3.49 


0.055 




MOON 


3.4 t of ""*Mo 


1 1027 


[21] 


2.78 


0.024 




CAMEO 


1 t of CdWOi crystals 


^ 1026 


[21] 


2.42 


0.085 




CUORE 


750 kg of TeOa 


^ 1027 


[197] 


2.95 


0.023 




XMASS 


10 t of liq. Xe 


3 1026 


[21] 


1.97 
1.67 


0.062 
0.073 
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Figure 43: Lower bound on neutrinoless double beta decay in the models of Refs. [132] (left) and [133] 
(right). In first model the lower bound depends on the Majorana phase (pi, while in the second one has 
a strong dependence of the lower bound on the value of the atmospheric angle. 

to an effective mass parameter M^e ^ 0.03 eV, as illustrated in the right panel of the same figure. 

Complementary information on the absolute scale of neutrino mass is expected to come from future 
beta decays searches such as in the KATRIN experiment [202]. Altogether beta and PPou sensitivities 
should be confronted with information coming from cosmology [68, 69]. 
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9 Beyond neutrino oscillations 



How robust is the current oscillation picture of neutrino data? For example, how well do we know the 
astrophysics of the Sun, neutrino propagation and/or neutrino detection to be confident that there are 
no loop holes in our understanding of neutrinos? 

9.1 Solar magnetic fields 

Here we give two examples of how solar magnetic fields can affect the neutrino oscillation interpretation 
of solar neutrino data. 

9.1.1 Radiative zone magnetic fields and density noise 

The Sun can harbor magnetic fields in its radiative zone as well as in its convective zone. It has been 
shown [218] that magnetic fields deep within the solar radiative zone can produce density fluctuations 
that could affect solar neutrino fluxes in an important way [219]. However it has been shown that 
consistency with KamLAND reactor data restores robustness of the determination of neutrino oscillation 
parameters [217]. 

9.1.2 Convective zone magnetic fields and spin fiavour precession 

Convective zone magnetic fields can lead to spin flavour precession if neutrinos have non-zero transition 
magnetic moments [28] that may affect neutrino propagation both in vacuo and in matter [220, 221]. A 
global analysis of spin-flavour precession (SFP) solutions to the solar neutrino problem is characterized 
by three effective parameters: Amg^L = Am^, the neutrino mixing angle ^sql = and the magnetic fleld 
parameter IjlB± [222, 223]. For n = 10^^^ Bohr magneton, and an optimum self-consistent magneto- 
hydrodynamics magnetic fleld proflle with maximum strength ~ 80 KGauss in the convective 
zone [224] one flnds an excellent description of the solar neutrino data in terms of spin flavor precession. 
However, after combining with data from the KamLAND reactor experiment, one flnds that such 
solutions to the solar neutrino problem are ruled out, as they can not account for the deflcit and 
spectral distortion observed at KamLAND. 

9.1.3 Probing for neutrino magnetic moments 

Even though by themselves Majorana neutrino transition moments can not provide an acceptable ex- 
planation of the solar phis KamLAND data, they can still be present at a sub-leading level and be 
studied using existing solar and reactor neutrino data. They would contribute to the neutrino-electron 
scattering cross section and hence alter the signal observed in Super-Kamiokande and at Borexino. In 
ref. [225] constraints were placed on the neutrino transition moments by using the solar neutrino data. 
It was found that all transition elements can be bounded at the same time. Furthermore, improved 
reactor data play a complementary role to the solar neutrino data in further improving the sensitivity. 
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which is currently at 2 x 10~^°/iB at the 90% CL. for the combined solar + reactor data. The upcoming 
Borexino experiment will improve the bounds from today's data by roughly one order of magnitude, 
thanks to the lower energy. 

9.2 Non-standard neutrino interactions 

It is hard to find a neutrino mass generation scheme that does not bring in dimension-6 non-standard 
neutrino interaction (NSI) terms. Such sub- weak strength sGp operators are illustrated in Fig. 44. They 
can be of two types: flavour- changing (FC) and non-universal (NU) and may arise in several ways. One 
is from the non-trivial structure of charged and neutral current weak interactions characterizing the 
broad class of seesaw-type models [7]. If present, these interactions would in general bring additional 



sources of CP violation. While the expected magnitude of the non-standard interactions is rather model 
dependent, it may well fall in the range that can be tested in the precision neutrino oscillation studies 
discussed in Sec. 6. The first manifestation non-standard interactions can have is through the violation 
of unitarity in the lepton mixing matrix [7]. The non- unitary piece of the lepton mixing matrix can 
be sizable in inverse seesaw-type models [43] and hence can be phenomenologically important [80, 81, 
82]. With neutrino physics entering the precision age it becomes an important challenge to scrutinize 
how good is the unitary approximation of the lepton mixing matrix in future experiments, given its 
theoretical fragility. 

Relatively sizable NSI strengths may also be induced through the exchange of new scalar bosons, 
as present in models with radiatively induced neutrino masses. These typically contain new relatively 
light scalar states such as Higgs bosons [50], scalar leptoquarks, etc. NSI terms may also arise from 
renormalization group evolution in supersymmetric unified models [71]. 

Non-standard physics may affect neutrino production and detection cross sections, as well as prop- 
agation properties. In their presence, the Hamiltonian describing neutrino propagation has, in addition 
to the standard oscillation part, another term i^NSb 



ward scattering amplitude for the FC process t'/i + / — > + / and ^/^GpNfe' represents the difference 




Figure 44: Flavour-changing effective operator for non-standard neutrino interaction. 




(101) 



Here +(— ) holds for neutrinos (anti-neutrinos) and e and e' parametrize the NSI: y^GpNfe is the for- 
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between 1^^ + f and i^r + f elastic forward scattering. Here Nf is the number density of the fermion 
/ along the neutrino path. A remarkable fact is that, in the presence of non-standard interactions 
resonant neutrino conversions can take place even in the absence of neutrino masses [203]. 

Solar neutrino oscillations 

Although non-standard interaction effects may provide an alternative way to account for the current 
solar neutrino data [204], consistency with KamLAND reactor neutrino data requires them to be sub- 
leading [113]. However the oscillation interpretation of solar neutrino data is still "fragile" against the 
presence of non-standard interactions in the e — r sector, opening another degenerate solution in the 
second octant of the solar mixing angle [205]. Indeed, as seen in Fig. 45, in the presence of NSI, there 
appear new solar neutrino oscillation solutions in addition to the standard one. The two degenerate 
solutions are denoted LMA-I or "normal" and LMA-D "degenerate or dark-side" solution. Although 




sin e^o^ sin e^o^ 

Figure 45: New solar neutrino oscillation solutions in the presence of NSI [205]. 



some ways of inducing the new LMA-D solution, for example via NSI between neutrinos and down-type- 
quarks-only, are already in conflict with the combination of current atmospheric data and accelerator 
neutrino oscillation data of the CHARM experiment, the general form of these solutions is not yet ruled 
out. Moreover, further precision KamLAND reactor measurements will not resolve the ambiguity in the 
determination of the solar neutrino mixing angle, since they are expected to constrain mainly ArngQ^- 

However, as can be seen from Fig. 46, the "normal" LMA-I solutions have quite different predicted 
neutrino survival probabilities for low energy neutrinos (left panel) and boron neutrinos (right panel). 
Thus future precision low energy solar neutrino experiments will help in lifting the degeneracy between 
the "normal" and the new degenerate solution. 

On the other hand, even if very small, non-standard interaction effects may play an important role in 
supernova astrophysics, opening the possibility of new resonances that could take place in the internal 
neutron-rich regime [206, 207, 208, 209]. 
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Figure 46: Probing for NSI in future low energy solar neutrino experiments [205]. 
Atmospheric neutrino oscillations 

In contrast to the solar case, it has been shown [210] that, within the 2-neutrino approximation, 
the determination of atmospheric neutrino parameters Am^^^j and sin^ 6'atm has been shown to be prac- 
tically unaffected by the presence of NSI on down- type quarks {f = d). There are however loop-holes 
in the three neutrino case [211]. Let us also mention here that future neutrino factories will have good 
potential for probing non-standard neutrino in the mu-tau channel [212]. 



NSI-oscillation confusion theorem 

For the coming generation of experiments we note that even a small residual non-standard interaction 
of neutrinos in the e — r channel characterized by a parameter eer can have dramatic consequences for 
the prospects of probing neutrino oscillations at neutrino factories. For example, it has been shown [213] 
that the presence of such NSI leads to a drastic loss in sensitivity in the ^13 determination at a neutrino 
factory. It is therefore important to improve the sensitivities on NSI, another window of opportunity 
for neutrino physics in the precision age. 



10 Summary 

We have reviewed the basic mechanisms to generate neutrino mass, analysing the corresponding struc- 
ture of the lepton mixing matrix. After an overview of the status of neutrino oscillation parameters 
as determined from current data, we have focused on the next neutrino oscillation experiments using 
accelerator neutrinos, and discussed the prospects for probing the strength of CP violation in two of 
them, T2K and NOz/A. Intermediate term experiments using wide band beams and capable of probing 
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also the second Oscillation Maximum, such as T2KK, have also been discussed [182, 184, 185, 183]. 

Farther in the future are experiments with novel neutrino beams, such as neutrino factories and 
beta beams. One of the advantages of the neutrino beam coming from a muon storage ring is that the 
energy spectrum is accurately known, in contrast with conventional muon neutrino beam from pion 
decay, where the energy spectrum is not very well determined. It was suggested in [214] that pure Ve 
or Pg beams can be obtained by using (3 unstable isotopes. The idea is to produce first a huge number 
of (5 unstable ions and then accelerate them in a storage ring to some reference energy, and let them to 
decay in the straight section of a storage ring in order to get an intense or Pg beam. The advantage 
of this method is that the energy spectrum is precisely known (as in the case of neutrino beam from 
muon decay) and the there is no contamination of the other neutrino species unlike the case of neutrino 
beam from muon decay. None of these possibilities were considered here, but they are discussed in a 
recent review by Geer and Zisman in this same series [215] for neutrino factories, and e.g., in Ref. [216] 
for the case of beta beam. Another extensive review is being prepared for the International scoping 
study of a future Neutrino Factory and super-beam facility. We have also briefly commented on the 
possibility of probing CP violation effects induced by Major ana phases in neutrinoless double beta decay. 
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